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ABSTRACT 


The  behavior  of  both  unbounded  and  bounded  compressible  plane 
mixing  layers  with  respect  to  two-  and  three-dimensional,  spatially  growing  wave 
disturbances  is  investigated* using  linear  stability  analysis.  The  mixing  layer  is 
formed  by  two  parallel  streams  with  different  gases  and  the  flow  is  assumed  to  be 
inviscid  and  non-reacting. 


For  unbounded  mixing  layers,  the  effects  of  the  free-stream  Mach  num¬ 
ber,  velocity  ratio,  temperature  ratio,  gas  constant  (molecular  weight)  ratio  and  the 
ratios  of  specific  heats  on  the  linear  spatial  instability  characteristics  of  a  mixing 
layer  are  determined.  A  nearly  universal  dependence  of  the  normalized  maximum 
amplification  rate  on  the  convective  Mach  number  is  found  for  two-dimensional 
spatially  growing  disturbances.  The  effects  of  the  mean  flow  profiles  on  the  insta¬ 
bility  behavior  of  the  mixing  layers  are  also  studied.  It  is  shown  that  decreasing 
the  thickness  of  the  total  temperature  profile  relative  to  tne  mean  velocity  profile, 
or  adding  a  wake  component  in  the  mean  velocity  profile  can  make  the  normalized 
amplification  rate  decrease  slower  as  the  convective  Mach  number  increases  for  both 
subsonic  and  supersonic  convective  Mach  numbers. 


For  an  unbounded  mixing  layer  with  subsonic  convective  Mach  num¬ 
bers,  there  is  only  one  unstable  mode  propagating  with  a  phase  velocity  C*,^  approx? 
imately  equal  to  the  isentropically  estimated  convective  velocity  of  the  large  scale 
structures  u*.  As  the  convective  Mach  number  approaches  or  exceeds  unity,  there 


are  always  two  unstable  spatial  modes. ^One  is  with  a  phase  velocity  C*m  <  u*  (slow 
mode)  and  the  other  is  with  a  phase  velocity  C*m  >  u*  (fast  mode).  For  the  low 
supersonic  convective  Mach  numbers,  the  fast  mode  is  more  unstable  than  the  slow 
mode  when  the  heavy  gas  is  on  the  low  speed  side  and  the  slow  mode  is  dominant 
when  the  heavy  gas  is  on  the  high  speed  side. 

The  effect  of  parallel  flow  guide  walls  on  a  spatially  growing  mixing 
layer  is  also  investigated.  It  is  shown  that,  in  this  case,  if  the  convective  Mach 
number  exceeds  a  critical  value  of  approximately  unity,  there  are  many  supersonic 
unstable  modes.  The  maximum  amplification  rates  of  mixing  layers  approach  an 
asymptotic  value  and  this  maximum  amplification  rate  increases  to  a  maximum 
value  and  decreases  again  as  the  distance  between  the  walls  decreases.  For  a  mixing 
layer  inside  parallel  flow  guide  walls,  the  growth  rate  of  three-dimensional  modes 
is  larger  than  the  corresponding  two-dimensional  mode  at  high  convective  Mach 
numbers.  But  the  growth  rate  of  two-dimensional  supersonic  instability  waves  has 
a  larger  value  than  their  three-dimensional  counterparts  for  a  mixing  layer  inside  a 
rectangular  duct  (Tam  &  Hu  [1988],  [1989]).  Contour  plots  of  the  pressure  pertur¬ 
bation  fields  for  both  unbounded  and  bounded  mixing  layers  indicate  that  there  are 
waves  propagating  outward  from  the  mixing  layer  along  the  Mach  angle,  and  that 
the  walls  provide  a  feedback  mechanism  between  the  growing  mixing  layer  and  this 
compression/expansion  wave  system.  The  bounded  mixing  layers  axe  more  unstable 
than  the  corresponding  free  mixing  layers  for  supersonic  convective  Mach  numbers. 
The  streaklines  of  the  flow  confirm  that  the  spreading  rate  of  the  mixing  layer  is 
unusually  small  for  supersonic  disturbances. 
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CHAPTER  1 


INTRODUCTION 


1.1  Historical  Background 


The  aims  of  the  study  of  hydrodynamic  stability  were  to  find  how 
and  when  laminar  flows  break  down  and  undergo  transition  to  turbulence  or  some 
other  laminar  flow.  The  basic  problems  of  hydrodynamic  stability  were  studied  and 
formulated  by  Helmholtz  [1868],  Kelvin  [1871],  Rayleigh  [1880]  and  Reynolds  [1883] 
at  the  end  of  the  nineteenth  century.  The  hydrodynamic  equations  were  linearized 
by  assuming  that  the  disturbances  of  the  flows  were  infinitesimally  small.  The 
method  of  normal  modes  was  used.  Vivid  descriptions  of  the  problems  were  given 
by  Reynolds  through  his  series  experiments  on  the  instability  of  flow  in  a  pipe. 


The  first  analytic  study  of  two-dimensional  hydrodynamic  stability 
was  made  by  Helmholtz  [1868]  and  Kelvin  [1871].  They  considered  the  basic  flow  of 
incompressible,  inviscid  fluids  in  two  horizontal  parallel  infinite  streams  of  different 
velocities  and  densities,  one  stream  above  the  other,  and  solved  the  problem  of  the 
instability  of  a  wavy  disturbance  over  the  surface  of  discontinuity.  It  is  now  called 
the  Kelvin-Helmholtz  instability.  Later,  Rayleigh  [1913]  studied  the  instability  of 
incompressible,  inviscid  parallel  flows  with  continuous  velocity  distributions  across 
the  flows.  An  important  result  of  his  studies  is  the  finding  that  the  existence  of  an 
inflection  point  in  the  basic  velocity  profile  is  a  necessary  condition  for  instability. 
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This  result  was  further  enhanced  by  the  work  of  Tollmien  [1935].  He  pointed  out 
that  this  “inflection-point  criterion”  is  not  only  necessary  but  also  sufficient  for 
flows  with  a  symmetrical  or  a  boundary-layer  type  velocity  distribution.  Further 
development,  clarification  and  physical  interpretation  of  these  general  results  were 
given  by  Lin  [1944,  1955]. 

As  a  result  of  the  increasing  importance  of  the  phenomena  of  gas 
flow  at  high  speeds,  linear  stability  investigations  have  been  extended  to  include 
the  effect  of  compressibility.  The  first  attempt  to  develop  a  compressible  stability 
theory  was  made  by  Kuchemann  [1938].  He  studied  the  stability  of  a  piecewise 
linear  velocity  profile  representing  a  boundary  layer  in  a  compressible  fluid.  Later 
on,  the  stability  of  tangential  discontinuities  in  a  compressible  fluid  was  treated  by 
Landau  [1944].  He  showed  that  tangential  discontinuities  axe  stable  with  respect  to 
infinitesimal  disturbances. 

Important  theoretical  investigations  of  the  stability  of  the  compress¬ 
ible  laminar  free  and  bounded  mixing  layers  to  infinitesimal  wave  disturbances  were 
made  by  Lees  &  Lin  [1946].  Their  study  was  in  the  form  of  an  extension  of  the 
principles  and  techniques  already  formulated  for  the  study  of  the  stability  of  incom¬ 
pressible  laminar  flows  to  compressible  laminar  flows.  They  gave  a  general  stability 
criterion  in  terms  of  the  gradient  of  the  product  of  density  and  vorticity,  analogous 
to  the  Rayleigh-Tollmien  criterion  for  the  case  of  an  incompressible  fluid.  The  dis¬ 
turbances  outside  the  mixing  layers  were  classified  by  Lees  &  Lin  as  “subsonic,” 
“sonic”  and  “supersonic”  disturbances,  according  to  whether  the  phase  velocity  of 
the  wavy  disturbance  in  the  direction  of  the  free  stream  relative  to  that  of  the  fluid 
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is  less  than,  equal  to,  or  greater  than  local  velocity  of  the  sound.  Also,  a  detailed 
consideration  to  an  inviscid  theory  of  compressible  mixing  layers  was  given  and  the 
possibility  that  the  disturbances  in  a  compressible  flow  can  produce  a  turbulent 
transition  or  lead  to  the  formation  of  a  shock  was  also  pointed  out  by  Lees  &:  Lin. 


The  instability  of  inviscid  compressible  free  mixing  layers  with  respect 
to  two-  and  three-dimensional,  temporally  growing  disturbances  was  considered  by 
Lessen,  Fox  &  Zien  [1965,  1966]  for  both  subsonic  and  supersonic  disturbances. 
Under  the  assumption  that  the  flow  was  iso-energetic,  h*  =  h*  +u*2/2  =  constant, 
they  found  that  the  flow  is  unstable  with  respect  to  supersonic  disturbances,  al¬ 
though  the  amplification  rate  is  smaller  than  for  subsonic  disturbances  and  the 
increase  of  the  angle  between  the  disturbance  wave  number  vector  and  the  principle 
flow  direction  tends  to  increase  the  instability.  With  spatially  growing  disturbances, 
Gropengiesser  [1970]  studied  this  instability  problem  using  the  Crocco-Busemann 
relation,  h*  =  h*  +  u*2 / 2  =  fn(u*),  for  the  mean  temperature  profile,  which  would 
be  discussed  in  detail  later,  and  using  the  compressible  laminar  boundary  layer  ve¬ 
locity  profile  for  the  mean  velocity  distribution  across  the  mixing  layers.  He  carried 
out  the  inviscid  instability  calculations  at  various  free  stream  Mach  numbers  and 
temperature  ratios.  He  found  the  existence  of  a  second  unstable  mode  and  showed 
that  a  change  in  the  character  of  the  disturbance  may  involve  the  amplification  of 
a  second  mode.  He  also  noted  the  high  growth  rate  of  three-dimensional  waves 
at  high  Mach  numbers.  In  order  to  simplify  the  stability  problems,  Blumen  and 
Drazin  et  al.  [1970,  1975,  1977]  studied  the  stability  problem  by  assuming  that  the 
thermodynamic  state  of  a  compressible,  inviscid,  free  mixing  layer  is  constant.  They 
showed  that  two-dimensional  disturbances  are  unstable  at  all  values  of  the  Mach 
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number  and  that  there  exists  a  second  unstable  supersonic  mode. 


A  critical  review  of  the  experimental  data  for  developed  compressible 
free  turbulent  mixing  layers  was  given  by  Birch  &  Eggers  [1973].  They  indicated 
that  although  some  inconsistency  exists  among  the  various  experimented  investiga¬ 
tors,  there  is  a  definite  trend  of  decreasing  growth  rate  with  increasing  free  stream 
Mach  number.  Experimental  evidence  of  large  scale  structures  in  two-stream  super¬ 
sonic  turbulent  mixing  layers  were  given  by  Chinzei  et  al.  [1986],  by  Papamoschou 
&  Roshko  [1986,  1988]  and  Papamoschou  [1986,  1989].  Papamoschou  &  Roshko 
suggested  the  convective  Mach  number  Mc  as  the  appropriate  parameter  scaling 
the  effect  of  compressibility.  This  is  defined,  for  each  stream,  as: 


Mcl  =  -1- . 


Mc2  = 


u;  —  u  2 
a2 


(1.1) 


where  uj,  uj  and  °i>  °2  3X6  the  free  stream  velocities  and  speeds  of  sound.  The 
quantity  u*  is  the  convective  velocity  of  the  large  scale  structures.  Coles  [1981] 
pointed  out  that  stagnation  points  must  exist  in  these  large  scale  structures.  The 
convective  velocity  u*  was  estimated  by  Papamoschou  &:  Roshko  assuming  that  the 
dynamic  pressures  match  at  stagnation  points  in  the  flow  (Coles  [1985],  Dimotakis 
[1986]).  For  compressible  isentropic  flow  with  equal  static  pressures  in  the  two  free 
streams,  e.g.,  Papamoschou  &  Roshko, 


(1.2) 


where  71,  72  are  the  ratios  of  the  specific  heats  of  the  two  streams.  If  71  and  72  are 
equal,  it*  is  given  by 


u 


* 

C 


a^u*  4-  aju2 
a\  +  a\ 


0 


(1.3) 
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which,  for  equal  static  free  stream  pressures,  reduces  to  the  incompressible  ex¬ 
pression  (Dimotakis  [1986]).  They  suggested  that  the  normalized  growth  rate  of 
a  compressible  mixing  layer  (unity  for  incompressible  flow)  might  be  expressible 
as  a  universal  function  of  the  convective  Mach  number  Mc\  which  is  valid  over  a 
wide  range  of  velocity  and  temperature  ratios  of  a  mixing  layer.  They  also  indi¬ 
cated  that  the  normalized  growth  rate  decreased  significantly  with  increasing  Mc\ , 
for  subsonic  Mc  1 ,  and  reached  an  asymptotic  value  for  supersonic  convective  Mach 
numbers.  The  same  qualitative  behavior  had  been  found  by  Bogdanoff  [1983]  in  his 
analysis  of  several  previous  experimental  investigations  of  supersonic  mixing  layers. 


Using  linear  spatial  instability  theory,  Ragab  &  Wu  [1988,  1989]  stud¬ 
ied  the  influence  of  velocity  ratio  on  the  stability  characteristics  of  compressible  free 
mixing  layer.  They  also  investigated  the  effect  of  the  convective  Mach  number  on 
the  growth  rates  of  mixing  layers.  Their  results  indicate  that  the  convective  Mach 
number  is  a  parameter  which  correlates  compressibility  effects  on  the  spreading 
rate  of  mixing  layers.  Jackson  &  Grosch  [1988,  1989]  presented  their  results  of  a 
study  of  the  inviscid  spatial  instability  behavior  of  compressible  free  mixing  layers 
with  one  stream  moving  and  the  other  stream  stationary.  They  showed  that  if  the 
Mach  number  of  the  moving  stream  exceeds  a  critical  value,  there  axe  always  two 
groups  of  unstable  waves.  One  of  these  groups  is  feist  with  normalized  phase  speeds 
greater  than  1/2,  and  the  other  is  slow  with  phase  speeds  less  than  1/2.  They 
also  showed  that  three-dimensional  modes  have  the  same  general  behavior  as  two- 
dimensional  modes,  but  with  higher  growth  rates  over  some  range  of  propagation 
direction.  Recently,  Sandham  &  Reynolds  [1989a,  1989b]  investigated  compressible 
mixing  layers  using  both  linear  theory  and  direct  simulation.  They  showed  that  lin- 
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ear  theory  can  be  very  useful  in  understanding  the  physics  of  free  mixing  layers  and 
the  growth  rates  of  the  developed  plane  mixing  layers  and  that  three-dimensional 
modes  are  dominant  in  an  unbounded  high-speed  mixing  layer  above  a  convective 
Mach  number  of  0.6. 


For  a  compressible  boundary -layer,  the  study  of  linear  instability  the¬ 
ory  carried  out  by  Mack  [1969, 1984]  indicated  that  for  supersonic  disturbances  there  • 

are  two  kinds  of  neutral  waves,  inflectional  and  noninflectional,  and  that  there  are 
always  unstable  modes  in  boundary-layers.  Mack  also  noted  the  high  amplification 

rate  of  three-dimensional  waves  at  high  Mach  number.  % 


For  a  mixing  layer  inside  a  rectangular  channel,  Tam  &:  Hu  [1988] 
showed  that  the  coupling  between  the  motion  of  the  mixing  layer  and  the  acous¬ 
tic  modes  of  the  channel  produces  new  instability  waves  for  the  spatially  growing 
mixing  layers.  Greenough  et  al.  [1989]  considered  the  effects  of  walls  on  a  confined 
compressible,  temporally  growing  mixing  layer.  They  showed  that  there  are  two 
general  types  of  instabilities:  confined  Kelvin- Helmholtz  modes  and  supersonic  wall 
modes.  A  class  of  highly  amplified  supersonic  disturbances  are  found  by  Macaraeg 
&  Streett  [1989]  for  high-speed,  temporally  growing  bounded  shear  flows  at  high 
values  of  the  streamwise  wave  number. 


1.2  Present  Research  Objective 


# 


i  ne  present 
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instability  behavior  of  both  unbounded  and  bounded  compressible  plane  mixing 
layers  with  respect  to  two-  and  three-dimensional,  spatially  growing  wave  distur¬ 
bances.  The  purpose  of  this  study  is  to  understand  the  effect  of  compressibility  on 
spatially  growing  plane  mixing  layers  and  also  the  effect  of  the  mean  flow  profiles 
on  the  instability  behavior  of  the  mixing  layers.  The  mixing  layer  is  formed  by  two 
parallel  streams  with  different  gases  and  the  flow  is  assumed  to  be  inviscid  and  non 
reacting. 

The  convective  velocity  is  estimated  by  the  phase  velocity  of  the  most 
unstable  mode  C*m,  since  we  only  consider  the  amplified  disturbances  (a,  <  0). 
Therefore,  the  convective  Mach  number  for  each  stream  is  defined  here  as: 


For  these  amplified  disturbances,  we  think  the  definition  given  by  Eq.  1.4  is  appro¬ 
priate  since  we  are  tracking  the  most  unstable  mode.  Mack  [1975]  used  the  phase 
velocity  of  the  neutral  mode  c*N  as  the  convective  velocity  in  his  stability  analysis 
of  the  existence  of  many  supersonic  neutral  waves  of  supersonic  boundary  layer. 
Sandham  &  Reynolds  [1989b]  proposed  that  the  large-scale  structures  found  in  the 
mixing  layer  are  associated  with  neutral  instability  modes.  Therefore,  they  used 
the  phase  velocity  of  the  neutral  mode  CN  to  estimate  the  convective  velocity  of  the 
large-scale  structures  in  the  mixing  flow. 

For  a  compressible  mixing  layer  comprised  of  different  gases,  the  in¬ 
stability  calculations  were  carried  out  for  inviscid  wave  disturbances,  t.e.,  the  effect 
of  viscosity  was  neglected  in  the  disturbances,  and  the  mean  flow  profiles  were 
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assumed  to  vary  smoothly  across  the  mixing  layer.  The  disturbances  were  also  as¬ 
sumed  to  have  small  amplitudes.  Therefore,  the  linearized  equations  could  be  used. 
The  main  flow  was  assumed  to  be  quasi-parallel,  i.e.,  it  possessed  only  a  velocity 
component  in  the  direction  of  the  flow.  These  assumptions  have  been  justified  by 
previous  investigators  ( e.g .  Lin  [1955]).  The  assumption  of  parallel  flow  provided 
an  even  better  approximation  in  the  case  of  compressible  mixing  layers,  since  the 
growth  rate  of  the  mixing  layer  became  smaller  with  increasing  free  stream  Mach 
number  M\. 


The  linearized  equations  for  inviscid  disturbances  were  derived  and 
the  eigenvalue  problem  was  formulated  in  Chapter  2.  The  range  of  unstable  frequen¬ 
cies  and  wave  numbers  was  numerically  calculated  using  a  Runge-Kutta  method 
combined  with  a  shooting  technique.  The  mean  flow  distributions  were  either  as¬ 
sumed  or  derived  from  the  boundary  layer  equations  for  compressible  mixing  flow. 
The  details  of  these  mean  flow  profiles  were  discussed  in  Chapter  3. 


For  compressible,  unbounded  mixing  layers,  the  influence  of  free 
stream  Mach  number,  velocity  ratio,  temperature  ratio  and  gas  constant  (molecular 
weight)  ratio,  and  ratios  of  the  specific  heats  on  the  linear  spatial  instability  char¬ 
acteristics  of  a  mixing  layer  were  investigated.  The  effect  of  the  convective  Mach 
number  on  the  growth  rate  of  the  mixing  layers  was  studied.  Numerical  calculations 
of  linear  spatial  instability  characteristics  of  the  mixing  layers  with  different  combi¬ 
nations  of  the  free  stream  Mach  number,  velocity  ratio,  temperature  ratio  and  gas 
constant  ratio  were  performed  in  which  the  free  stream  Mach  number  ranged  from 
0  to  11,  the  velocity  ratio  varied  from  0.25  to  0.75,  the  temperature  ratio  varied 
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from  0.2  to  5  and  the  gas  constant  ratio  varied  from  0.1  to  5.  In  order  to  assess 
the  sensitivity  of  the  details  of  the  mean  flow  profiles  on  the  instability  behavior  of 
the  compressible  mixing  layers,  an  initial  mean  velocity  profile  with  a  wake  compo¬ 
nent  and  a  total  temperature  profile  with  a  different  thickness  relative  to  the  mean 
velocity  profile  were  also  considered. 


Finally,  for  a  mixing  layer  inside  parallel  flow  guide  walls,  the  effects 
of  walls  and  the  distance  between  the  walls  on  a  spatially  growing  mixing  layer  were 
investigated.  For  a  given  distance  between  the  parallel  flow  guide  walls,  the  linear 
spatial  instability  characteristics  of  mixing  layers  were  calculated  for  different  free 
stream  Mach  numbers.  Contour  plots  of  the  pressure  perturbation  fields  indicate 
that  there  are  waves  propagating  outward  from  the  mixing  layer  along  the  Mach 
angle,  and  that  reflections  of  the  compression/expansion  waves  caused  by  walls 
provide  a  feedback  mechanism  between  the  growing  mixing  layer  structures  and  the 
wave  system  for  supersonic  convective  Mach  numbers.  Also,  the  flow  patterns  of 
the  mixing  layers  were  obtained  by  calculating  the  streaklines  of  the  mixing  layers. 
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CHAPTER  2 


MATHEMATICAL  FORMULATION 


The  objectives  of  this  chapter  axe  to  derive  the  linearized  disturbance 
equations  for  inviscid,  compressible  plane  mixing  layers  and  to  formulate  the  eigen¬ 
value  problem  of  the  instability  of  the  mixing  layers.  The  mixing  layer  is  formed 
by  two  parallel  streams  with  different  perfect  gases  and  is  subjected  to  a  two- 
dimensional  spatially  growing  wave  disturbance. 


2.1  Equations  for  Small  Disturbances 


Consider  a  two-dimensional  compressible  mixing  layer  formed  by  two 
parallel  streams.  Assume  that  the  fluids  on  each  stream  of  the  mixing  layer  are 
inviscid  and  non  reacting.  In  the  absence  of  external  forces,  the  general  equations 
of  a  two-dimensional,  inviscid,  compressible  flow  of  a  perfect  gas  can  be  written  in 
Cartesian  co-ordinates  as  follows: 


(a)  Equation  of  continuity: 


Dp *  ,  0JdvT_  d£\ 

Dt*  \3j*  dy* ) 


=  0 


(b)  Equations  of  motion: 

du *  tdu*  mdu *  1  dp* 

-  +  u - 1-  v  -  - - — 

dt*  dx*  dy*  p*  dx* 


(2.1) 


(2.2) 
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dv *  ,dv*  ,dv *  1  dp * 

dt*  dx*  dy *  p*  dy * 


(2.3) 


(c)  Equation  of  energy: 


/ dT *  ,ar* 


dT' 


du*  dv* 


+  «•  a?  +  •' V )  =  -p' W ' +  W )  {2A) 

(d)  Equation  of  state: 

p*  =  p*R*T *  .  (2.5) 

All  quantities  in  these  equations  Eire  dimensioned,  eis  defined  in  the  List  of  Symbols. 


# 


Since  the  compressible  mixing  layers  we  considered  are  formed  by  two 
different  ideed  gases,  according  to  the  Gibbs-Dalton  law,  the  pressure  in  a  mixture 
of  different  gases  is  equal  to  the  sum  of  the  partial  pressures  of  the  different  gases, 
i.e., 

P*  =  P*R*T*  =  rf  +  rf  =  p\R\T*  +  p*R*T*  ,  (2.6) 

and  the  corresponding  expression  for  enthalpy  can  be  written  as 

h*  =  XI K  +  x*h*2  =  (xic;,  +  X2C;2)T* ,  (2.7) 

where  xi  is  the  mass  fraction  of  species  1  of  the  gas  mixture.  The  mass  fraction, 
X2  >  of  species  2  is  equal  to  1  —  xij  since  only  two  species  exist  in  the  mixing  layer. 
Thus  Eqs.  2.6  and  2.7  yield 

ir  =  xii*r  +  (l-Xi)i*2*,  (2.8) 

c;  =  xic;,  +a  -  xi)c;2.  (2.9) 

For  perfect  gases,  a  similar  expression  can  be  obtained  for  the  specific  heat  at 
const  suit  volume  c*,  i.e., 

=  X1K1  +  (1  -  X1K2  • 


(2.10) 
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Since  the  mixing  layer  is  formed  by  two  different  gases,  Eqs.  2.1  to  2.5  do  not  form 
a  complete  set.  They  have  to  be  supplemented  by  the  species  transport  equation, 
which,  in  the  absence  of  molecular  diffusion,  is  given  by 


dXi  ,  ,dXi  ,  . dxi  n 
dr  +U  dx*  +V  dy •  0  ‘ 


(2.11) 


We  consider  a  time-independent  basic  flow  and  a  small  amplitude, 
time-dependent  disturbance.  For  a  parallel  flow,  the  basic  steady  state  is  of  the 


form 


Q *  =  <?V)eIt 


(2.12) 


where  ex  denotes  the  unit  vector  in  the  x*  direction.  The  general  quantity  of  the 
flow  can  be  expressed  as 


Q*(x*,y*,n  =  QV)  +  QV>yV*), 


(2.13) 


where  Q  ( y )  is  a  mean  flow  profile,  which  will  be  discussed  in  Chapter  3,  and  Q'  is  a 
small  amplitude  disturbance.  By  substituting  Eq.  2.13  for  each  of  the  variables  into 
Eqs.  2.1  to  2.5  and  Eq.  2.11  and  neglecting  terms  quadratic  in  the  small  disturbance, 
we  obtain  the  linearized  disturbance  equations.  With  upper  stream  quantities  as  the 
reference  and  the  local  layer  thickness  6*(x*)  as  the  length  scale,  the  dimensionless 
variables  are  defined  as  follows: 

t  =  x  =  x*/ 6\  y  =  y*/6\ 

U  =  U*/u  1*,  V  =  V*  fu\, 

p  =  p-/Pl\  T  =  T*/T*,  R  =  *•//?,*, 
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p  —  P  fp\  i  Cp  —  Cp  /Cpi  ,  Cu  —  Ct,  /c„i  . 

The  non  dimensional  parameters  of  the  problem  are  the  upper  stream  Mach  number 
Mi  =  u*/a*,  where  a*  is  the  speed  of  sound  of  the  upper  stream,  and  the  ratio 
of  specific  heats  of  the  upper  stream  71  =  cjh/cJh-  The  linearized  dimensionless 
disturbance  equations  are  arrived  at  : 


Continuity: 


Momentum: 


dp'  _  dp'  ,dp  _  ( du'  dv'  \ 

+"3r  +  !,3;+p(a7  +  ^)  =  0 

1  dp' 


( du '  du'  .du 

ibr+i^  +  t,VI  =  - 


dx 


dy 


Jdv'  _dv'\ 


\  at 


+  udl)  =  - 


7xMi2  dx 
1  dp' 


liMy2  dy 


(2.14) 

(2.15) 

(2.16) 


Energy: 


State: 


(dv  dr  ,dT\  ,  ,jdu'  dv'\  /01„, 


Species  transport: 


1  =  pRT,  p'  =  ^  +  =■  +  £ 
P  R  T 


d\\  d\\  3y/ 

+  tH^-  +  =  0. 

at  dx  dy 


(2.18) 


(2.19) 


In  these  equations,  u,  p  axe  dimensionless  mean  flow  quantities  and  u',  ...,  p' 
are  the  corresponding  dimensionless  disturbances.  Eqs.  2.14  to  2.19  form  a  sys¬ 
tem  of  dimensionless  linear  partial  differential  equations  for  the  disturbances.  The 
coefficients  are  known  functions  given  by  the  steady  mean  flow  profiles.  After  lin¬ 
earization,  the  dimensionless  expressions  of  Eqs.  2.8  to  2.10  become: 

~  #2  .  R?\_  /,  i?2\. 


(2.20) 


1 
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(2.21) 

(2.22) 


2.2  Normal  Mode  Solution  for  a  Spatially  Growing  Mixing  Layer 


We  consider  here  the  disturbance  to  be  a  two-dimensional  wave  prop¬ 
agating  in  the  x-direction.  The  disturbance  quantities  in  a  suitably  defined  dimen¬ 
sionless  form  can  be  expressed  as: 

u',v'  =  {/(y),a^(y)}  exp  [  ia(x  -  ct)  ]  ;  (2.23) 

p',T',p'  =  {r(y),9(y),n(y)}  exp  [  m(x  -  ct)  ]  ;  (2.24) 

Xi  =  Tj(y)  exp  [  ia(x  -  ct)  ] ;  (2.25) 

where,  for  a  spatially  growing  layer,  a  is  a  dimensionless  complex  wave  number 
and  c  is  a  dimensionless  complex  wave  velocity.  Each  component  travels  with  the 
dimensionless  phase  speed  Cp  =  u>/ar  (u>  =  ac)  and  grows  or  dies  away  in  x  like 
exp(— Oix).  Therefore  the  disturbances  are  amplified,  neutral,  or  damped,  according 
to  whether  a*  <  0,  «<  =  0,  or  <*{  >  0,  respectively.  It  is  to  be  understood  that  the 
real  part  of  the  disturbances  represents  the  physical  quantity  in  each  case. 
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2.3  Differential  Equations  for  the  Complex  Amplitude  Functions 


Using  normal  mode  solutions,  we  can  separate  the  variables  and  re¬ 
duce  the  linearized  disturbance  equations  from  the  partial  differential  equations  to 
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ordinary  differential  equations.  By  substituting  Eqs.  2.23  to  2.25  into  Eqs.  2.14 
to  2.19,  we  obtain  a  system  of  ordinary  differential  equations  for  determining  the 
amplitude  functions  f(y),  p(y ),  &(y),  7r (y)  and  rj(y).  These  are: 


Continuity: 


Momentum: 


i(u  —  c)r  +  4>~f  +  p(if  4-  <t>)  =  0  , 
ay 


,du 

,4>- T  +  ~  c)f 

[  dy 


pa2i(u  —  c)<f>  —  — 


ITT 


7i  Mi2  ' 


Energy: 


State: 


pcv 


i(u  —  c)8  +  <f> 


dT 

dy 


7i  M?  ' 

=  —  (7i  -  1)  (*/  +  <t>)  , 


Species  Transport: 


ir  r  6  1  ,  R.2  \ 

p  ~  ~p  +  T  +  R^1  ~  R^V  ’ 


i(u  -  c)r i  +  <t> ^  =  0  , 


(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 


where  dots  denote  d/dy.  These  equations  for  the  amplitude  functions  can  be  reduced 
to  the  following  second-order  differential  equation  for  the  pressure  disturbance. 


7T  — 


2  5  _  (T  +  R\ h  _a2  L_  ~  c)2cv/R 

u-c  \T  Rj\  \  rtL/r  +  (7i  _  i) 


►7T  =  0  .  (2.32) 


All  the  other  variables  can  be  expressed  in  terms  of  it  and  if  sis  follows: 


-/  = 


7!  Ml2 p(u  -  c) 


IT  + 


TT  du 


a2(u  —  c)2  dy  J 


(2.33) 
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7i  Mi2a2p(u  —  c) 


(2.34) 


2.4  Boundary  Conditions 


Two  types  of  boundary  conditions  will  be  considered.  One  is  for 
unbounded  mixing  layers  and  the  other  is  for  bounded  mixing  layers.  These  will  be 
discussed  below. 


2.4.1  Unbounded  Mixing  Layer 


The  general  configuration  of  a  free  mixing  layer  and  the  variables 
involved  are  summarized  in  Fig.  2.1.  The  boundary  conditions  of  the  eigenfunction 
7r (y)  are  found  from  Eq.  2.32.  With  y  — ►  ±oo,  the  mean  quantities  u(y),  T(y),  R(y) 
and  c„(y)  are  constants  and  their  corresponding  derivatives  are  zero.  In  that  limit, 


Eq.  2.32  becomes 


7 r  —  A±27t  =  0 


(2.35) 


A+2  =  q2[1  -  M2(  1  -  c)2]  =  A+  =  A+r  +  tA+i 


a2  (l  -  7iMi2(ur  -  c)2cvr/Rr  )  _  A  _  A 
^  RrTr  Cvr/Rr  +  (jl  -  1)  ' 

>  a 


(2.36) 


—  A_r  +  iA-i  ,  (2.37) 
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and  the  solution  for  large  |y|  are  of  the  form 

7T  =  A±  exp  (  ±A ±y  )  , 
where  is  a  complex  constant. 


(2.38) 


Since  we  are  only  considering  the  case  of  amplified  spatially  growing 
disturbances  (cti  <  0),  the  boundary  conditions  for  both  subsonic  and  supersonic 
disturbances  require  that  7r(y)  must  be  bounded  as  y  — ♦  ±oo.  To  define  A±  uniquely, 
we  set  the  real  part  of  A±=A±r  >  0  to  obtain 


y  =  v+ oo  -*•  +00 
y  =  y- oo  -oo 


tt  =  A+  exp  (  -A +y  )  , 
n  =  A-  exp  (  +A_y  )  , 


(2.39) 

(2.40) 


where 

A±  =  A±r  +  i  A±, 


n  l1/2  fi  11/2 

=  -dA±|  +  A±.)  +tsign(A±i)  -(|A±|-A±r)  .(2.41) 


2.4.2  Bounded  Mixing  Layer 

The  general  configuration  of  a  mixing  layer  inside  parallel  walls  is 
shown  in  Fig.  2.2.  6  is  the  dimensionless  mixing  layer  thickness  defined  in  a 

boundary-layer  sense,  i.e., 

<f  and  u2---"(-*)  <e,  (2.42) 


where  e  is  taken  as  0.002.  Note  that  2d  is  the  dimensionless  distance  between  the 
walls. 
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In  the  regions  of  6  <  y  <  d  and  —  d  <  y  <  — 6 ,  the  flow  is  effectively 
uniform,  i.e.,  the  mean  flow  quantities,  such  as  velocity,  temperature,  etc.  are 
constants,  and  Eq.  2.32  reduces  to  Eq.  2.35.  In  order  to  satisfy  the  boundary 
conditions  at  the  walls,  the  vertical  velocities  of  the  disturbances  should  vanish  on 
the  solid  boundaries  and  we  must  have  tt (d)  =  -k(— d)  =  0  .  The  analytic  solutions 
in  Regions  1  and  2  then  are 

7r  =  A  cosh  [  A+(<f  —  y)],  tt  =  B  cosh  [  A_(d  +  y)  ]  ,  (2.43) 

where  A+  and  A_  are  already  given  by  Eqs.  2.36  and  2.37.  A  and  B  are  complex 
constants  to  be  determined  by  matching  to  the  inner  flow  region. 


2.5  Formulation  and  Numerical  Treatment  of  the  Eigenvalue  Problem 


The  eigenvalue  problem  of  the  instability  of  mixing  layers  with  spa¬ 
tially  growing  disturbances  is  defined  as  follows.  We  take  a  read  value  of  the  distur¬ 
bance  frequency  w  and  seek  the  complex  eigenvalue  a  =  ar  +  ia,.  The  dispersion 
relation  can  be  written  in  the  form 

F(a,u>)  =  0.  (2.44) 

For  a  given  real  disturbance  frequency  u>,  the  complex  eigenvalue  a  is  to  be  deter¬ 
mined  such  that  the  eigenfunction  n (y)  of  Eq.  2.32  satisfies  the  boundary  conditions. 
This  gives  a  well-defined  eigenvalue  problem  for  the  amplified  disturbances  a,  <  0. 
A  Runge-Kutta  method  combined  with  a  shooting  technique  was  used  to  solve  this 
eigenvalue  problem. 
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For  an  unbounded  mixing  layer,  Eq.  2.32  is  integrated  numerically 
for  a  fixed  value  ui  starting  from  the  lower  boundary  y-oo-,  where  the  boundary 
condition  is  given  by  Eq.  2.40,  to  the  upper  boundary  y+oo.  This  yields  7r(y+00) 
and  tt (y+0O).  By  matching  the  upper  boundary  condition  given  by  Eq.  2.39,  the 
sum  S(ar,0!j)  =  7r(y+0O)  +  A+7r(y+0O)  is  evaluated  and  improved  values  of  a  are 
calculated  from  the  approximated  zeros  of  S(ar,ati)  by  linear  interpolation.  This 
procedure  is  repeated  until  1 5|  is  sufficiently  small. 


For  a  mixing  layer  inside  parallel  flow  guide  walls,  we  use  the  analytic 
0  solution  in  Region  2  of  Fig.  2.2  as  a  starting  solution  and  integrate  Eq.  2.32  numer¬ 

ically  from  y  =  —6  to  y  =  6.  The  correct  a  is  obtained,  for  a  given  u>,  by  matching 
the  numerical  solution  to  the  analytic  solution  at  y  =  6. 


The  computations,  for  both  free  and  bounded  mixing  layers,  were 
performed  using  the  Runge-Kutta-Fehberg  method  with  automatic  estimation  of 
local  error  and  step  size  adjustment.  An  error  control  of  10-5  is  used  for  the 
integrations,  iterating  until  the  eigenvalue  converges  to  10~6. 
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CHAPTER  3 


MEAN  FLOW  QUANTITIES 


In  a  finite  thickness  mixing  layer,  all  mean  quantities,  such  as  veloc¬ 
ity,  density,  temperature,  gas  constant  (molecular  weight),  concentration  and  the 
ratio  of  the  specific  heats,  vary  gradually  across  the  mixing  layer.  We  consider  the 
two-dimensional  motion  of  a  stream  of  fluid  with  the  velocity  u*,  density  /?*,  tem¬ 
perature  T*,  gas  constant  /?*,  ratio  of  the  specific  heats  7*,  viscosity  /x*  and  thermal 
conductivity  k *,  over  a  parallel  stream  described  by  the  corresponding  parameters 
with  subscript  2.  The  main  flow  is  considered  to  be  quasi-parallel  and  the  fluids  in 
the  two  parallel  streams  are  assumed  non-reacting. 


A  two-dimensional  compressible  steady  mixing  flow,  comprised  of  dif¬ 
ferent  perfect  gases,  can  be  modelled  by  the  compressible  boundary  layer  equations 
for  gas  mixture.  We  take  the  axis  of  x*  to  be  in  the  direction  of  motion  of  the  free 
streams,  and  the  axis  of  y*  to  be  in  the  transverse  direction.  The  velocity  compo¬ 
nents  in  these  directions  are  u*  and  v*.  The  origin  is  taken  as  the  point  at  which 
the  two  fluids  first  come  in  contact  (see  Fig.  2.1).  For  a  steady  parallel  mixing  flow, 
the  disturbance  term  Q'  in  Eq.  2.13  is  zero,  therefore,  Eq.  2.13  becomes 


=  (3.D 

With  zero  pressure  gradient  throughout  the  entire  flow  region,  the  global  boundary - 
layer  equations  for  compressible  steady  mixing  flow  are 
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Continuity  equation: 

0(p*O  d(p*v*) 

dx *  dy* 

Momentum  equation: 


(3.2) 


(3.3) 


Energy  equation: 


'■(•■S  *  •■£) 


Equation  of  state  for  ideal  gases: 


p*  =  p*R*T* 


(3.5) 


Species  transport  equation: 


.  .dXi  ,  .  ,5Xi  _  d  f  ,5xi  \ 

P  U  dx*  +  P  V  dy*  “  dy*  V°  P  dy*  )  ’ 


(3-6) 


where  p*  is  the  viscosity,  Pr  =  p*c*/k*  is  the  Prandtl  number,  Sc  =  p*/p*D*  is 
the  Schmidt  number,  k *  is  the  thermal  conductivity,  and  D*  is  the  binary  diffusion 
coefficient. 


3.1  Mean  Velocity  Distribution 

The  numerical  calculation  of  the  mean  velocity  distribution  for  an 
incompressible  mixing  layer  with  u^/ui  ^  0,  comprised  of  the  same  gases  was  first 
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provided  by  Lock  [1951].  By  introducing  a  stream  function  ip*,  which  satisfies  the 
continuity  equation  of  compressible  steady  flow  Eq.  3.2  identically,  i.e., 


.  sr 

W' 


.  .  dr 

p  v  = 


dx*  ’ 


(3.7) 


and  a  transformation  proposed  by  L.  Howarth  [1948], 

r'  =  /f  ■ 


(3.8) 


the  calculation  of  the  mean  velocity  distribution  for  compressible  steady  mixing 
flows  is  converted  into  the  corresponding  calculation  for  incompressible  flows,  if  we 
assume  that  the  viscosity  p*  of  the  gas  is  approximately  a  linear  function  of  the 
temperature  T*  (Gropengiesser  [1970]).  It  was  found  by  Gropengiesser  [1970]  that 
the  dimensionless  velocity  distribution  u(y)  for  compressible  steady  mixing  flows 
can  be  well  approximated  by  the  generalized  hyperbolic  tangent  profile  with  three 
free  constants, 


u(y)  =  1  ~  ~  (!  ~  tanh  nx  (y  -  yx))  ] 


KJ 


(3.9) 


with  ki  =  0.307257,  «2  =  3.695640,  y*  =  2.127137.  The  dimensionless  variable  y  is 
equal  to  Y* /6*(x*),  where  <5*(x*)  is  the  local  layer  thickness.  Therefore,  the  main 
flow  is  considered  to  be  quasi-parallel. 


In  order  to  simplify  the  problem,  for  the  mixing  flows  considered  in 
the  current  numerical  calculations,  we  assume  that  the  dimensionless  mean  velocity 
profile  is  described  by  a  hyperbolic  tangent  profile  of  the  form 


«(y)  =  C(y)  +  «r[i  -  C(y)] , 


(3.10) 


where  ur  =  U2/U1  is  the  velocity  ratio  across  the  mixing  layer,  and 

C(y)  =  |[1  +  tanh(y)]  .  (3.11) 

See  mean  velocity  profiles  u(y)  with  different  values  of  ur  in  Fig.  3.1.  It  has 
been  known  from  previous  investigators  that  the  percentage  difference  in  ampli¬ 
fication  rate  of  the  case  with  a  hyperbolic  tangent  profile  relative  to  the  case  with  a 
boundary-layer  profile  is  small  if  the  value  of  the  velocity  ratio  ur  is  not  too  small. 


It  is  known  that  there  is  a  wake  component  in  the  initial  evolution  of 
the  mixing  layer  (e.g.  Lang  [1985]).  For  the  incompressible  mixing  layer,  Kooches- 
fahani  &  Frieler  [1987]  investigated  the  linear  spatial  instability  behavior  of  a  mean 
velocity  distribution  with  a  wake  component.  They  showed  the  existence  of  two 
unstable  modes,  which  they  called  shear  layer  mode  and  wake  mode.  In  order  to 
investigate  the  effect  of  the  wake  component  on  the  instability  behavior  of  the  com¬ 
pressible  mixing  layer,  we  assumed  that  the  mean  dimensionless  velocity  profile  is 
composed  of  the  usual  hyperbolic  tangent  profile  plus  a  wake  component  (owing  to 
the  splitter  plate)  represented  by  a  Gaussian  distribution  of  the  form 

«(y)  =  C(y)  +  Mi  -  C(y)]  -  ^  exp(-y2) ,  (3.12) 

where  w  is  the  normalized  wake  deficit,  as  was  assumed  in  Koochesfahani  &  Frieler 
[1987].  Profiles  with  different  values  of  w  are  plotted  in  Fig.  3.2. 


3.2  Mean  Temperature  Distribution 


Two  different  kinds  of  temperature  profiles  have  been  considered.  For 
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the  flow  comprised  of  same  gases,  according  to  Crocco  [1932]  and  Busemann  [1935], 
it  is  possible  to  satisfy  the  energy  equation  identically  by  assuming  that  the  tem¬ 
perature  T*  is  dependent  only  upon  the  velocity  component  u*  taken  parallel  to 
the  direction  of  the  flow.  Following  the  Crocco  and  Busemann’s  idea,  we  assume 
that  the  total  enthalpy  h *  =  h*  +  u*2/2  depends  only  on  the  velocity  component  u* 
for  the  mixing  flows  comprised  of  different  gases,  i.e.  h*  =  fn(u*).  By  substituting 
the  expressions  h *  =  A*(u*)  into  the  energy  equation  3.4  and  by  assuming  that  the 
Prandtl  number  Pr  =  p*  c*/k*  and  Schimdt  number  Sc  =  p*/p*D *  are  equal  to 


unity,  i.e., 

Pr  =  1 
Sc  =  1, 

and  using  the  momentum  Eq.  3.3,  Eq.  3.4  can  be  reduced  to 


(3.13) 


dPh •  _  cP(h*  +  m*2/2) 
du *2  du*2 


(3.14) 


Thus  the  mean  static  temperature  profile  can  be  represented  by 

J?(V*)T*(W*)  =  +  c,  «•(, ,•)  +  a  (3.15) 

where  cx  and  ci  are  constants  determined  by  the  boundary  conditions  on  the  mean 
temperature  profile  and  we  also  assumed  that  c*  ^  fn{T*)  within  the  range  of  the 
temperatures  in  the  flow  and  c*  is  only  a  function  of  the  fraction  of  species  of  the 
gas  mixture.  The  dimensionless  mean  temperature  profiles  T(y)  with  a  fixed  A/j, 
or  a  fixed  temperature  ratio  T2/T1,  are  shown  in  Figs.  3.3  and  3.4,  respectively. 
Equations  3.10  and  3.15  yield  the  dimensionless  toteil  temperature  profile  Tt(y)  for 
the  mixing  flows  formed  by  seune  gases, 


r,(»)-((»)  +  T,r[  1  -  <(»)], 


(3.16) 
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where  Ttr  =  Tt2/T<i  is  the  ratio  of  the  total  temperature  across  the  mixing  layer. 
The  profile  given  by  Eq.  3.16  has  the  same  thickness  relative  to  the  dimensionless 
mean  velocity  profile  given  by  Eq.  3.10.  A  question,  then  remains  as  to  how  this 
relative  thickness  affects  the  instability  behavior  of  the  compressible  mixing  layer. 
To  consider  this  question,  £(y)  is  replaced  by  C(y/<7)  in  Eq.  3.16,  where  a  adjust  the 
thickness  of  the  total  temperature  profile  relative  to  the  velocity  profile.  The  dimen¬ 
sionless  total  temperature  profiles  as  a  function  of  the  total  temperature  interface 
thickness  a  are  shown  in  Fig.  3.5. 

An  alternate  dimensionless  mean  temperature  profile  considered  was 
obtained  by  assuming  that  the  profile  can  also  be  represented  by  a  hyperbolic  tan¬ 
gent  profile,  i.e., 

T(»)-C(»)  +  T,[l  -«»)],  (3.17) 

where  Tr  =  T2/Ti  is  the  ratio  of  the  temperature  across  the  mixing  layer.  Com¬ 
paring  Eq.  3.17  with  Eq.  3.15,  we  see  that  the  profile  given  by  Eq.  3.17  is  only  a 
function  of  the  temperature  ratio  Tr.  Therefore,  for  both  subsonic  and  supersonic 
flows,  we  have  the  same  dimensionless  mean  temperature  profiles  for  the  same  Tr 
(see  Fig.  3.6). 

3.3  Mean  Concentration  Distribution 

Since  the  mixing  layer  we  considered  may  be  formed  by  different  gases, 
the  mean  mass  fractions  and  *2  may  vary  gradually  across  the  mixing  layer, 
just  like  the  mean  velocity  and  temperature.  Following  the  way  we  specify  the 
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dimensionless  mean  velocity,  we  assume  that  the  mean  concentration  profile  \i(y) 

is  also  described  by  a  hyperbolic  tangent  profile  of  the  form 

Xi(y)  =  j[l  +  tanh(y)]  . 

(3.18) 

Since  Xi  4-  X2  =  1>  Eq.  3.18  yields 

A 

X2 (y)  =  \il  -  tanh (y)]  . 

(3.19) 

Substituting  Eq.  3.18  into  Eqs.  2.20  to  2.22,  the  dimensionless  mean  gas  constant 

and  specific  heat  profiles  can  be  written  as: 

*(y)  =  C(y)  +  *r[i  -  CO)] , 

(3.20) 

cP(y)  =  CO)  d*  cpr  [  i  -  C(y)] , 

(3.21) 

cv(y)  =  CO)  +  cvr[l  -  C(y)] , 

(3.22) 

where  C(y)  was  given  by  Eq.  3.11,  Rr,  and  cpr,  cvr  are  the  ratio  of  the  gas  constant 
and  the  ratios  of  the  specific  heats  at  constant  pressure  and  at  constant  volume 
respectively.  The  dimensionless  mean  profile  R(y )  is  shown  in  Fig.  3.7  for  different 
values  of  Rr. 
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CHAPTER  4 


THE  INSTABILITY  BEHAVIOR  OF  A  COMPRESSIBLE 


FREE  MIXING  LAYER 


In  this  chapter,  the  linear,  spatial  instability  behavior  of  a  compress¬ 
ible,  laminar,  mixing  layer,  comprised  of  different  gases,  is  investigated  with  respect 
to  two-  and  three-dimensional,  spatially  amplified,  wave  disturbances.  The  fluctu¬ 
ating  amplifying  flow  is  assumed  to  be  inviscid  and  non  reacting.  The  perturbations 
in  the  flow  are  taken  as  possessing  small  amplitudes.  All  flow  quantities  mentioned 
in  this  chapter  are  defined  as  dimensionless  quantities,  using  the  dimensionless  vari¬ 
ables  defined  in  Chapter  2.  The  effects  of  free  stream  Mach  number,  velocity  ratio, 
temperature  ratio,  gas  constant  (molecular  weight)  ratio  and  the  ratio  of  specific 
heats  on  the  linear  spatial  instability  characteristics  of  the  compressible  mixing  lay¬ 
ers  are  examined.  The  behavior  of  eigenfunctions  of  the  disturbed  mixing  layers  is 
discussed  for  both  subsonic  and  supersonic  disturbances.  The  instability  character¬ 
istics  of  compressible  mixing  layers  with  spatially  growing  oblique  wave  disturbances 
are  also  studied.  The  influence  of  adding  a  wake  component  to  the  mean  velocity 
profile  and  decreasing  the  thickness  of  the  total  temperature  profile  relative  to  the 
mean  velocity  profile  on  the  linear  spatial  instability  behavior  of  compressible  mix¬ 
ing  layers  is  also  considered. 
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4.1  Free  Stream  Conditions 

The  instability  characteristics  of  compressible  mixing  layers  can  be 
calculated  numerically  using  Eq.  2.32,  with  given  boundary  conditions  and  mean 
flow  conditions.  The  mean  flow  distributions  of  mixing  layers  were  already  described 
in  Chapter  3.  For  a  plane  mixing  layer  having  the  mean  velocity  profile  Eq.  3.10, 
the  mean  temperature  profile  Eq.  3.15  and  the  mean  concentration  profile  Eq.  3.18, 
the  instability  behavior  of  mixing  layers  with  two-dimensional,  spatially  growing 
disturbances  was  studied  in  the  range  of  free  stream  Mach  number  0  <  M\  <  11.0, 
velocity  ratio  0.25  <  u2/ui  <  0.75,  temperature  ratio  0.2  <  T2/T\  <  5.0,  gas 
constant  ratio  0.1  <  R2/R\  <  5.0  for  given  ratios  of  specific  heats  71,  72.  The  lists 
of  the  combinations  of  the  free  stream  conditions  described  above  are  presented  in 
Tables  4.1  -  4.3. 

The  mean  temperature  distribution  across  the  mixing  layer  is  also 
approximated  by  a  hyperbolic  tangent  profile  Eq.  3.17.  The  mean  flow  conditions 
for  this  hyperbolic  tangent  mean  temperature  profile  are  given  in  Table  4.4. 
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Table  4.1  The  free  stream  conditions  for  mixing  flows  formed  by  the  same  gases 


Table  4.2  The  free  stream  conditions  for  different  values  of  R\  and  R2 


Table  4.4  The  free  stream  conditions  for  hyperbolic  tangent  mean  temperature 
distribution 

Ri/R-i  =  1.0.  7i  =  1.4,  71  =1.4 


Case 

M, 

u2/ui 

t2/t1 

18 

0  —  4.5 

0.50 

0.67 

19 

0  -►  5.0 

0.50 

1.00 

20 

0  -►  5.5 

0.50 

2.00 

The  instability  characteristics  were  calculated  for  each  combina¬ 
tion  of  the  mean  flow  conditions,  yielding  the  most  unstable  eigenvalue  am 
( am  =  arm  +  i  Oim),  its  corresponding  frequency  u>m  and  phase  velocity 

Cpm  =  ^m/^rm  •  (4.1) 

The  convective  Mach  numbers  Me  1 ,  Mc 2 ,  corresponding  to  these  spatially  growing 
waves,  can  be  then  obtained  from  Eq.  1.4. 

4.2  Effect  of  the  Free  Stream  Mach  Number  Mi 

The  inviscid  solutions  of  the  instability  characteristics  of  compressible 
mixing  layers  were  calculated  for  various  free  stream  Mach  numbers  M\  for  a  fixed 
velocity  ratio  u2/ui,  temperature  ratio  T2/T1,  gas  constant  ratio  R2/R1  and  ratio 
of  the  specific  heats  72. 

For  a  free  mixing  layer  with  subsonic  convective  Mach  numbers,  there 
is  only  one  group  of  unstable  waves.  The  amplification  rate  —a,  is  shown  in 
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Figs.  4.1a,  4.2a  and  4.3a  for  Cases  2,  15  and  4  (see  Tables  4.1  and  4.2),  as  a  function 
of  the  amplified  disturbance  frequency  u>  and  the  free  stream  Mach  number  Mi .  The 
maximum  amplification  rate  — oim  attains  its  biggest  value  for  an  incompressible 
mixing  layer  (Mi  =  0)  and  falls  off  sharply  with  increasing  free  stream  Mach  num¬ 
ber  Mj .  The  band  of  amplified  frequencies  becomes  much  narrower  with  increasing 
free  stream  Mach  number  M\.  These  results  are  in  agreement  with  Gropengiesser’s 
[1970]  numerical  calculations,  who  considered  a  laminar  velocity  profile  with  the 
velocity  ratio  ui/ui  —  0.  It  should  also  be  noticed  in  Figs.  4.1a,  4.2a  and  4.3a  that 
for  a  constant  free  stream  Mach  number  M\ ,  the  amplification  rate  —a,  increases 
monotonically  with  increasing  amplified  disturbance  frequency  w,  up  to  a  maximum 
value,  and  then  drops  off  gradually  to  zero  as  the  neutral  disturbance  frequency  is 
reached.  The  corresponding  phase  velocities  Cp  =  uj/ar  of  these  instability  char¬ 
acteristics  are  given  in  Figs.  4.1b,  4.2b  and  4.3b  against  the  amplified  disturbance 
frequency  u  for  different  free  stream  Mach  numbers  M\.  The  phase  velocity  Cp 
decreases  with  increasing  amplified  disturbance  frequency  u,  down  to  a  minimum 
value,  then  increases  slightly  till  the  neutral  disturbance  frequency  is  reached.  For 
a  mixing  layer  with  subsonic  disturbances,  the  phase  velocity  of  the  most  unstable 
wave  CPm  is  shown  in  Fig.  4.4  as  a  function  of  the  free  stream  Mach  number  M\ . 
From  this  figure,  we  can  see  that  Cpm  is  independent  of  the  free  stream  Mach  num¬ 
ber  Mi .  It  is  approximately  equal  to  the  isentropically  estimated  convective  velocity 
uc  defined  in  Eq.  1.3,  which  depends  only  on  the  velocity  ratio  temperature 

ratio  T2/T1,  density  ratio  P2/P1  and  the  gas  constant  ratio  R2/R1  for  71  =  72  (see 
Fig.  4.4). 


As  the  free  stream  Mach  number  M\  approaches  or  exceeds  a  critical 


value  Micr,  the  mixing  flow  becomes  supersonic  relative  to  its  disturbances,  which 
means  the  convective  Mach  number  Mc\  or  Mc 2  is  greater  than  unity.  Jackson  Sc 
Grosch  [1988]  showed  that  for  a  single  stream  (u2/ui  =  0)  supersonic  free  mixing 
layer,  there  are  always  two  groups  of  unstable  waves.  One  is  slow  with  a  scaled 
phase  velocity  less  than  1/2,  and  the  other  is  fast  with  a  scaled  phase  velocity 
greater  than  1/2.  They  classified  these  two  unstable  waves  as  a  slow  supersonic 
mode  and  a  fast  supersonic  mode  respectively.  The  linear  instability  investigations 
performed  here  indicate  that  for  a  supersonic  free  mixing  layer  with  the  velocity 
ratio  U2/U1  ^  0,  there  are  also  two  unstable  waves.  One  is  with  the  phase  velocity 
Cpm  <  uc  (slow  mode)  and  the  other  is  with  the  phase  velocity  Cpm  >  uc  (fast 
mode),  where,  as  noted  above,  uc  is  the  isentropically  estimated  convective  velocity 
(Eq.  1.3).  The  convective  Mach  number  Mc  1  is  always  greater  than  unity  and  Mc2 
is  usually  smaller  than  unity  for  the  slow  supersonic  mode  and  vice  versa  for  the 
fast  supersonic  mode. 

The  instability  characteristics  of  supersonic  mixing  flows  were  calcu¬ 
lated  for  both  unstable  supersonic  waves.  The  amplification  rate  — is  plotted  in 
Figs.  4.5a,  4.6a  and  4.7a  for  Cases  2,  15  and  4  at  a  given  free  stream  Mach  number 
Mi  >  Mi  cr •  The  bands  of  the  amplified  disturbance  frequencies  ui  of  the  two 
supersonic  unstable  waves  overlap  for  some  range  of  the  frequencies.  Therefore,  two 
unstable  modes  can  exist  at  the  same  time  for  some  given  amplified  disturbance  fre¬ 
quencies.  The  corresponding  variation  of  the  phase  velocity  is  given  in  Figs.  4.5b, 
4.6b  and  4.7b.  The  phase  velocities  of  the  most  unstable  mode  of  these  two  super¬ 
sonic  waves  Cpm  are  shown  in  Figs.  4.8a,  4.9a  and  4.10a  as  a  function  of  the  free 
stream  Mach  number  M\  >  M\cr,  for  Cases  2,  15  and  4.  From  these  figures,  we 


understand  that  the  phase  velocity  of  the  most  unstable  mode  Cpm  decreases  for 
a  slow  mode  and  increases  for  a  fast  mode  with  increasing  free  stream  Mach  num¬ 
ber  M\.  Therefore,  the  convective  Mach  number  Mc\  or  MC2 ,  defined  by  Eq.  1.4, 
increases  as  the  free  stream  Mach  number  increases.  The  maximum  amplification 
rates  —ctim  of  both  supersonic  unstable  modes  are  shown  in  Figs.  4.8b,  4.9b  and 
4.10b  as  a  function  of  the  free  stream  Mach  number  M\  >  Micr.  For  the  low  super¬ 
sonic  disturbances,  the  fast  mode  (Mc 2  >  1)  is  more  unstable  than  the  slow  mode 
(Mcl  >  1)  when  the  heavy  gas  is  on  the  low  speed  side  (Case  2)  and  the  slow  mode 
(Mc  1  >  1)  is  dominant  when  the  heavy  gas  is  on  the  high  speed  side  (Case  4  and 
Case  15).  This  might  suggest  that  if  there  is  a  shock  wave  in  the  mixing  flow,  the 
shock  would  be  borne  by  the  heavy  gas  stream  for  the  low  supersonic  convective 
Mach  numbers.  If  the  convective  Mach  numbers  are  high  enough,  the  maximum 
amplification  rate  of  the  two  supersonic  modes  are  comparable.  Therefore,  the  flow 
might  support  shocks  in  both  streams  for  the  high  convective  Mach  numbers. 


For  an  incompressible  mixing  layer  with  a  Blasius  mean  velocity  pro¬ 
file,  Monkewitz  &  Huerre  [1982]  found  that  the  maximum  amplification  rate  —  a,m 
is  nearly  proportional  to  the  velocity  ratio  A,  defined  as 

A  =  (1  -u2/ui)/(l  +  u2/ui)  .  (4.2) 

Instead  of  using  the  Blasius  mean  velocity  profile,  the  hyperbolic  tangent  mean 
velocity  profile  Eq.  3.10  was  used  in  the  current  calculations.  Other  mean  flow 
profiles  of  the  mixing  layers  were  obtained  from  Eqs.  3.15  and  3.18.  Sandham  & 
Reynolds  [1989b]  showed  that  the  differences  in  the  velocity  profile  are  generally 
small  when  the  two  streams  have  equal  densities. 
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In  order  to  study  the  effect  of  the  free  stream  Mach  number  M\  on 
the  function  —  ajm(A),  the  maximum  amplification  rate  — o,m  was  computed  as  a 
function  of  the  velocity  ratio  A  for  the  density  ratio  p2/Pi  =  1.0,  temperature  ratio 
Ta/Tx  =  1.0  and  the  ratio  of  specific  heats  71  =  72  =  1.4,  at  different  free  stream 
Mach  numbers  My .  The  result  shown  in  Fig.  4.11  indicates  that  the  relation  between 
the  maximum  amplification  rate  —  aim  and  the  velocity  ratio  A  is  approximately 
linear  for  incompressible  mixing  flow  with  equal  density  and  temperature  in  both 
mixing  streams  and  that  this  linear  relation  becomes  nonlinear  as  the  free  stream 
Mach  number  M\  increases. 


4.3  Effect  of  the  Velocity  Ratio  U2/111 


It  is  known  from  turbulent  mixing  layer  experiments  that  the  mix¬ 
ing  layer  spreading  rates  are  significantly  smaller  for  supersonic  mixing  flows  than 
those  observed  in  subsonic  experiments.  Although  the  effect  of  compressibility  on 
the  growth  rate  is  the  major  cause  for  the  lower  spreading  rate,  which  will  be  dis¬ 
cussed  later  in  Section  4.7,  the  effect  of  other  flow  parameters,  such  as  the  velocity 
ratio  u2/«i  should  also  be  studied.  Using  linear  spatial  instability  analysis,  the  ef¬ 
fect  of  the  velocity  ratio  u2/tti  on  the  instability  characteristics  of  mixing  layers  was 
considered.  For  a  temperature  ratio  T2/7j  =  0.5,  gas  constant  ratio  R2/R1  =  1.0 
and  a  ratio  of  specific  heats  71  =  72  =  1.4,  the  instability  characteristics  of  the 
mixing  layers  with  different  velocity  ratios  u2/ui  (Cases  2,  6  and  10)  are  shown 
in  Figs  4.12a,  4.12b,  4.12c  and  4.12d  at  different  convective  Mach  numbers  A/cl 
or  Mc 2,  which  can  be  obtained  from  Eq.  1.4.  For  both  subsonic  and  supersonic 
convective  Mach  numbers,  the  band  of  unstable  frequencies  u /  becomes  wider  and 


the  maximum  amplification  rate  —  atm  becomes  smaller  as  the  velocity  ratio  u2/u\ 
increases.  The  phase  velocities  Cpm  of  these  unstable  waves  are  plotted  in  Fig.  4.13 
against  the  velocity  ratio  u^/ui.  Fig.  4.13  demonstrates  that  the  phase  velocity  Cpm 
is  approximately  a  linear  function  of  the  velocity  ratio  u2/ui,  at  a  given  convective 
Mach  number  Mcl,  or  Mc 2,  and  other  free  stream  parameters  T2/Ti,  R2/R\,  71 
and  72.  The  solid  line  in  Fig.  4.13  is  the  isentropic  estimate  of  the  convective  veloc¬ 
ity,  Eq.  1.3,  uc(u2/ui),  which  is  independent  of  the  convective  Mach  number  Mc\ 
or  Mc 2.  For  subsonic  convective  Mach  numbers,  the  slope  of  the  linear  function 
Cpm(u2/ui)  is  the  same  as  the  slope  of  the  linear  function  uc(u2/uj).  For  super¬ 
sonic  disturbances,  there  are  two  groups  of  unstable  waves.  Even  though  the  linear 
relationship  between  the  phase  velocity  Cpm  and  the  velocity  ratio  u2/ui  exists  for 
each  unstable  wave,  the  slope  of  the  function  Cpm{u2lui )  is  greater  than  that  of 
the  isentropic  uc{u2/ul)  for  the  slow  supersonic  mode  ( Cpm  <  uc)  and  is  smaller 
for  the  fast  supersonic  mode  ( Cpm  >  uc). 

4.4  Effect  of  the  Density  Ratio  pzl  P\  or  the  Temperature  Ratio  T2/7i 

Plane  turbulent  mixing  between  two  parallel  streams  of  different  den¬ 
sities  was  studied  experimentally  by  Brown  k  Roshko  [1974]  for  incompressible 
mixing  layers.  The  effect  of  density  on  turbulent  mixing  in  incompressible  flow  was 
investigated.  A  comparison  of  the  results  for  the  growth  of  the  vorticity  thickness 
in  incompressible  flow  against  density  ratio  with  the  results  of  Maydew  k  Reed 
[1963]  in  supersonic  flow  was  made  in  Fig.  15  of  Brown  k  Roshko  [1974].  In  the 
current  investigations  of  the  linear  instability  behavior  of  inviscid  mixing  flows,  lam¬ 
inar  mixing  of  the  two  parallel  streams,  comprised  of  the  same  gas  (R2/Ri  =  1.0. 


7i  =  -y2),  was  considered.  The  effect  of  density  on  laminar  mixing  in  both  subsonic 
and  supersonic  flows  was  studied.  Figure  4.14  shows  the  maximum  amplification 
rate  —aim  as  a  function  of  the  density  r«.tio  P\/ P2  for  different  free  stream  Mach 
numbers  M\  at  the  velocity  ratio  ui/u\  =  0.  As  the  density  ratio  p\/p2  increases, 
the  maximum  amplification  rate  —  ajm  decreases  monotonically  for  an  incompress¬ 
ible  mixing  flow  (A/j  =  0).  For  a  compressible  mixing  flow,  —  a,m  increases  up  to  a 
maximum  value,  then  decreases  gradually  with  increasing  density  ratio  pi/p2-  The 
value  of  the  density  ratio  p\/p2,  which  corresponds  to  the  maximum  value  of  the 
maximum  amplification  rate  — a,m,  shifts  to  larger  values  of  the  density  ratio  pi/p2 
as  the  free  stream  Mach  number  M\  increases.  The  experimental  results  of  the 
growth  in  both  incompressible  and  compressible  mixing  flows  by  Brown  &  Roshko 
and  by  Maydew  &  Reed  are  also  given  in  Fig.  4.14.  The  numerical  calculations, 
however,  were  based  on  linear  instability  theory  while  the  experimental  results  were 
obtained  from  the  experiments  of  turbulent  mixing  layers.  We  can  see  that  the 
trend  of  the  instability  calculations  is  comparable  to  that  of  the  experimental  data. 

It  was  shown  in  Section  4.2  that,  for  incompressible  mixing  layers 
with  equal  density  on  both  mixing  streams  and  a  hyperbolic  tangent  mean  velocity 
profile,  the  maximum  amplification  rate  — a^m  is  approximately  a  linear  function  of 
the  velocity  ratio  A  =  (1  —  u2/ui)/(l  +  tt2/iti).  For  incompressible  mixing  layers,  the 
effect  of  the  density  ratio  p2/pi  on  the  maximum  amplification  rate  —  a,m  is  shown 
in  Fig.  4.15.  For  a  density  ratio  p2/Pi  equal  to,  or  less  than,  unity,  the  maximum 
amplification  rate  —  a,m  is  still  approximately  a  linear  function  of  the  velocity  ratio 
A  and  —  ctjm  rises  with  increasing  density  ratio  P2I P\  for  0  <  A  <  1.  For  a  density 
ratio  P2I p\  greater  than  unity,  the  function  — a,m(A)  becomes  nonlinear  and  the 
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maximum  amplification  rate  —  aim  falls  with  increasing  density  ratio  P2I  P\  for  some 
values  of  the  velocity  ratio  A.  For  incompressible  mixing  flows,  Dimotakis  [1986, 
1989]  proposed  an  expression  for  the  spatially  growing  mixing  layer  growth  rate 
given  by 


—  (ur,pr\Mi 
X 


(l-ur)(l+//2)[  (l-^)/(l +p‘/2)  ) 

2(1  \  1  + 2.9(1  +  ur)/(l  —  ur)  ]’ 

where  the  coefficient  Cg  is  independent  of  the  velocity  ratio  ur  and/or  the  density 

ratio  pr .  A  comparison  of  the  current  numerical  results  for  pr  =  7,  1/7  with  the 

experimental  values  of  Brown  &  Roshko  [1974]  and  with  the  Dimotakis  [1986,  1989] 

prediction  for  the  spatially  growing  mixing  layer  at  the  value  of  Cg  =  0.16  is  made 

in  Fig.  4.16. 


It  was  shown  in  Section  4.3  that  the  phase  velocity  of  the  most  unsta¬ 
ble  mode  Cpm  of  the  mixing  layer  is  approximately  a  linear  function  of  the  velocity 
ratio  u2/ui  at  a  given  convective  Mach  number  Mcl  or  Mc2  and  other  free  stream 
parameters  T2/T\,  R2/ R\,  and  j2.  Since  the  temperature  ratio  T2/T\  is  one  of  the 
important  parameters  for  mixing  layers,  the  effect  of  changing  the  temperature  ratio 
T2/T\  on  the  phase  velocity  Cpm  was  examined.  For  a  mixing  layer,  comprised  of  the 
same  gas  (R2/R\  =  1.0,  71  =  72  =  1.4),  the  phase  velocity  Cpm  was  calculated  as  a 
function  of  the  velocity  ratio  u2/u!,  at  a  given  T2/T\,  and  Afcl  or  Mc2.  The  effect  of 
the  temperature  ratio  T2/T1  on  the  function  Cpm(u2/ui)  is  shown  in  Fig.  4.17a  for 
subsonic  disturbances.  The  slope  of  the  approximately  linear  function  Cpm(u2/u  1) 
decreases  as  the  temperature  ratio  T2/T\  increases.  Also,  it  should  be  noted  that 
the  phase  velocity  Cpm  increases  with  increasing  temperature  ratio  T2/T1,  at  a  fixed 
velocity  ratio  u2/ui-  For  a  mixing  layer  {R2/Rx  =  1.0,  7j  =  72  =  1.4)  with  su¬ 
personic  convective  Mach  numbers,  there  are  two  supersonic  unstable  modes:  the 


slow  supersonic  mode  with  the  phase  velocity  Cpm  <  uc,  and  the  fast  supersonic 
mode  with  the  phase  velocity  Cpm  >  uc.  The  functions  Cpm(u2/Ui)  for  the  two 
supersonic  unstable  modes  are  given  in  Fig.  4.17b  and  4.17c  respectively.  Compared 
with  the  subsonic  mode,  the  slope  of  the  approximately  linear  function  Cpm(u2/u\ ) 
is  greater  for  the  slow  supersonic  mode  Cpm  <  uc  and  is  smaller  for  the  fast  super¬ 
sonic  mode  Cpm  >  uc.  The  effect  on  the  function  Cpm(it2/ui)  of  each  supersonic 
mode,  however,  is  the  same  as  that  of  the  subsonic  mode  when  the  temperature 
ratio  T2/T1  is  changed. 

4.5  Effect  of  the  Ratios  of  Specific  Heats  7lt  72 

In  experiments  on  gas  phase  turbulent  mixing  layers,  a  layer  may  be 
formed  by  two  gases  with  different  ratios  of  specific  heats  7j  and  72.  The  effect  of 
the  ratios  of  specific  heats  71 ,  72  on  the  instability  characteristics  of  mixing  layers 
could  be  important.  Here  we  studied  this  effect  using  linear  instability  analysis. 

The  instability  characteristics  of  mixing  layers  were  calculated  for 
Cases  7,  16  and  17.  For  incompressible  mixing  layers,  the  instability  characteristics 
are  not  affected  by  the  difference  between  71  and  72  (see  Fig.  4.18a).  For  subsonic 
mixing  layers  with  Mc\  0,  the  band  of  the  amplified  disturbance  frequencies 
becomes  wider  and  the  maximum  amplification  rate  —  aim  increases  if  the  high¬ 
speed  stream  has  the  lower  value  of  the  ratio  of  specific  heats  (see  Fig.  4.18b). 
The  effect  of  the  ratios  of  specific  heats  71,  72  on  the  instability  behavior  of  the 
slow  supersonic  mode  Cpm  <  uc  is  the  same  as  for  the  subsonic  unstable  mode 
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(see  Fig.  4.19a).  On  the  other  hand,  the  band  of  amplified  disturbance  frequencies 
becomes  narrower  and  the  maximum  amplification  rate  — otm  becomes  smaller  for 
the  fast  supersonic  mode  Cpm  >  uc,  when  the  high-speed  stream  has  the  lower  value 
of  the  ratio  of  specific  heats  (see  Fig.  4.19b). 


4.6  Eigenfunction  Behavior 


For  a  given  eigenvalue  a,  the  eigenfunction  n (y)  for  the  pressure  dis¬ 
turbances,  corresponding  to  amplified  disturbances  (a,  <  0),  can  be  evaluated  by 
integrating  Eq.  2.32.  The  eigenfunctions  of  subsonic  and  supersonic  mixing  flows 
were  calculated  for  Case  2.  The  eigenfunction  ir (y)  of  the  most  amplified  eigenvalue 
am  is  shown  in  Fig.  4.20a  and  4.20b  for  subsonic  mixing  layers  with  free  stream 
Mach  numbers  Mi  — ♦  0  and  Mi  =  2.0,  respectively.  The  solid  lines  in  these  figures 
are  the  real  part  of  n(y)  and  the  dashed  lines  are  the  imaginary  part  of  7r(y).  As  can 
be  seen  from  the  equation,  both  reeil  and  imaginary  parts  of  the  eigenfunctions  n(y) 
decay  exponentially  away  from  the  mixing  layers  in  both  the  high-  and  low-speed 
streams.  For  supersonic  mixing  layers,  the  pressure  eigenfunction  n(y)  of  the  most 
amplified  eigenvalue  am  radiates  into  the  high-speed  stream  for  the  slow  supersonic 
mode  Cpm  <  uc  and  radiates  into  the  low-speed  stream  for  the  fast  supersonic  mode 
Cpm  >  uc.  But  for  these  two  supersonic  unstable  modes,  the  pressure  eigenfunc¬ 
tions  7r(y)  display  an  oscillatory  decay  in  the  stream  with  supersonic  disturbances 
and  an  exponentially  decay  in  the  stream  with  subsonic  disturbances  (see  Fig.  4.21a 
and  4.21b). 
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4.7  The  Universal  Dependence 

The  paper  by  Bogdanoff  [1983]  and  the  experimental  results  on  tur¬ 
bulent  mixing  layers  by  Chinzei  et  al.  [1986]  and  by  Papamoschou  &  Roshko  [1986, 
1988]  suggest  that  the  isentropically  estimated  convective  Mach  number  (defined  by 
Eq.  1.1)  is  the  appropriate  parameter  for  scaling  the  effect  of  compressibility.  Having 
defined  the  convective  Mach  number  of  the  turbulent  mixing  layers,  Papamoschou 
&  Roshko  argued  that  the  expression  of  the  growth  rate 

S'  =  fn(— ,  Mj, M2, 71,72)  (4.4) 

ui  p  1 

can  be  reduced  into  the  simpler  expression 

6'  =  fn(— ,^,7l,72,Mci)  ,  (4.5) 

u  1  Pi 

if  compressibility  effects  can  be  expressed  in  terms  of  a  single  parameter  Mcl .  After 
uncoupling  the  effect  of  the  convective  Mach  number  A/Cl,  they  suggested  that  the 
growth  rate  of  a  compressible  mixing  layer  might  be  generally  expressed  as 

~~F(Afci),  (4.6) 

°0 

where  S'Q  is  the  growth  rate  of  the  corresponding  incompressible  mixing  layer  and 
F  is  a  nearly  universal  function,  over  a  wide  range  of  velocity  ratios,  density  ratios, 
and  ratios  of  specific  heats. 


In  the  numerical  investigations  of  the  linear  instability  of  an  inviscid. 
compressible  plane  mixing  layer,  we  considered  a  two-dimensional,  spatially  growing 


wave  disturbance  of  the  form 
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where  a  is  a  complex  wave  number  and  c  is  a  complex  wave  velocity.  The  distur¬ 
bance  propagates  as  where  ui  is  a  real  disturbance  frequency.  The 

linear  instability  characteristics  of  mixing  layers  with  mean  flow  profiles  described 
in  Eqs.  3.10,  3.15  and  3.18  were  calculated  for  the  mean  flow  conditions  given  in 
Tables  4.1  -  4.3.  The  phase  velocity  of  the  most  unstable  mode  Cpm  was  obtained 
from  Eq.  4.1.  This  yields  the  convective  Mach  numbers  Mc\  and  Mc 2  from  Eq.  1.4. 


It  was  shown  in  Section  4.2  that  there  are  always  two  unstable  waves 
for  supersonic  mixing  layers.  One  is  with  the  phase  velocity  Cpm  <  uc  (slow  mode), 
and  the  other  is  with  the  phase  velocity  Cpm  >  uc  (fast  mode).  Since  the  slow 
supersonic  mode  Cpm  <  uc  is  more  unstable  than  the  fast  supersonic  mode  Cpm  > 
uc  for  most  cases  given  in  Tables  4.1  -  4.3,  only  the  slow  mode  was  considered  for 
supersonic  convective  Mach  numbers.  Results  shown  in  Figs.  4.22a  -  4.22h,  which 
were  obtained  from  different  mean  flow  conditions  listed  in  Tables  4.1  -  4.3,  indicate 
that  if  the  most  unstable  eigenvalue  for  a  compressible  mixing  layer  is  normalized 
by  its  value  corresponding  to  an  incompressible  mixing  layer  (at  the  same  velocity 
ratio,  density  ratio  and  the  ratio  of  the  specific  heats),  the  ratio  is  approximated  as 
a  function  of  the  convective  Mach  number  only  for  the  free  stream  conditions  given 
in  Tables  4.1  -  4.3.,  t.e.. 


S' (Mc  1) 
S'(Mci  =  0) 


max<  -a,(u2/ui,p2/pi,7i,72,Mci)  ^ 

- p - L-  ~  F(Md). 

max|-Qi(u2/ui,/?2/pi,7i,72,Mci  =  0)  [• 


(4.8) 


•  The  solid  line  estimate  of  6'(Mcl)/6'(Mc  1  =  0)  in  Figs.  4.22a  -  4.22h  was  computed 

using  all  the  data  corresponding  to  the  seventeen  different  mean  flow  conditions 
listed  in  Tables  4.1  -  4.3,  and  least  squares  fitting  the  normalized  maximum  ampli- 

•  fication  rate  vs.  the  convective  Mach  number  Mc\,  for  the  range  of  Mc\  from  0  to 
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2  with  a  function  of  the  form  (see  Fig.  4.23) 


Pci) 
6'(Mcl  =  0) 


1  +  Po 


e  — (pjMji  +  P3Atfi +P4Mj,)  _  j 


1  +  XPl 


(4.9) 


where 

po  =  0.758  ,  Pi  =  1.902  ,  P2  =  0.523  , 

Pi  =  —2.196  ,  P4  =  3.710  . 

Note  that  6'(MC i  — >  oo)/6'(Mci  =  0)  — ►  0,  i.e.,  the  normalized  maximum  am¬ 
plification  rate  of  unbounded  mixing  layers  decreases  continuously  with  increasing 
convective  Mach  number  in  the  supersonic  region.  Note  also  that  these  results 
suggest  that  ( dF/dMci)Mcl=o  =  0,  as  might  have  been  argued  a  priori.  Finally, 
we  note  that  the  quality  of  the  fit,  assuming  the  pi  =  2  (exactly),  which  yielded 
slightly  different  values  for  the  remainder  of  the  fit  parameters,  was  essentially  the 
same.  The  results,  shown  in  Figs.  4.22a  -  4.22h,  also  suggest  that  the  normalized 
maximum  amplification  rate  decreases  significantly  with  increasing  Mcl  for  subsonic 
convective  Mach  numbers.  In  the  region  1  <  Mc i  <  2,  this  normalized  amplifica¬ 
tion  rate  decreases  slowly,  but  continuously  as  the  convective  Mach  number  Mc\  is 
increased. 


In  the  second  set  of  calculations,  the  mean  temperature  profile  was 
assumed  to  be  a  hyperbolic  tangent  profile  and  was  obtained  from  Eq.  3.17.  The 
normalized  amplification  rate  was  calculated  for  different  values  of  the  temperature 
ratio  (Cases  18,  19  and  20)  as  a  function  of  the  convective  Mach  number  Mc\ 
from  0  to  about  2.0,  for  a  velocity  ratio  ui/ui  =  0.5  and  a  gas  constant  ratio 
Ri/R\  =  1.0.  The  results,  shown  in  Fig.  4.24,  substantiate  the  convective  Mach 
number  as  the  relevant  compressibility  parameter  and  also  display  good  agreement 


with  the  plot  of  6'(Mci)/6'(Mci  =  0)  obtained  from  Eq.  4.9,  even  though  these  two 
mean  temperature  profiles  are  very  different  at  supersonic  convective  Mach  numbers 
(see  Figs.  3.4  and  3.6). 


A  comparison  of  the  estimate  of  6'(Mc\)/8'{Mci  =  0)  with  Ragab  & 
Wu’s  numerical  data  and  with  Papamoschou  Roshko’s  experimental  data  is  made 
in  Fig.  4.25.  The  result  of  our  linear  instability  calculations  is  in  close  agreement 
with  Ragab  &  Wu’s  result,  which  was  also  based  on  linear  instability  calculations, 
but  using  a  laminar  velocity  profile  at  the  velocity  ratio  U2/U1  —  0.  Also,  Ragab 
&  Wu  used  the  isentropically  estimated  convective  velocity  uc  to  approximate  the 
convective  velocity  of  the  spatially  growing  wave  disturbances  instead  of  the  phase 
velocity  of  the  most  unstable  mode  Cpm,  as  was  done  here.  The  difference  between 
Cpm  and  uc,  though  not  small  in  the  region  of  supersonic  convective  Mach  numbers, 
does  not  affect  the  results,  since  the  normalized  amplification  rates  are  very  small  in 
this  region.  According  to  Papamoschou  &  Roshko’s  experimental  data,  the  growth 
rate  of  the  mixing  layer  approaches  an  asymptotic  value  as  the  convective  Mach 
number  becomes  supersonic.  As  opposed  to  the  experimental  findings,  however,  the 
growth  rate  of  the  unbounded  shear  layers,  based  on  linear  instability  calculations, 
decreases  continuously  in  the  supersonic  convective  Mach  number  region. 


4.8  Oblique  Wave  Disturbances 


The  numerical  calculations  described  above  were  carried  out  by  as¬ 
suming  that  the  disturbances  in  the  mixing  layers  are  two-dimensional  and  spa- 
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tially  growing.  In  the  work  of  Gropengiesser  [1970],  Ragab  &  Wu  [1988],  Jackson 
&  Grosch  [1988]  and  Sandham  Sz  Reynolds  [1989],  it  was  found  that  a  mixing  layer 
with  three-dimensional  wave  disturbances  has  the  same  general  behavior  as  that 
with  two-dimensional  wave  disturbance  but  with  higher  growth  rates  over  some 
range  of  the  propagation  direction  for  the  high-speed  mixing  layers. 


If  a  mixing  layer  is  subjected  to  a  small,  three-dimensional,  spatially 
growing  wave  disturbance  of  the  form 

q'(x,y,z,t)=q(y)e'(az  +  0*-“t\  (4.10) 

where  q  is  the  disturbance  amplitude,  a  and  j3  are  the  complex  wave  numbers  in 
the  downstream  x  and  cross-stream  2  directions  respectively,  and  u>  is  the  real  dis¬ 
turbance  frequency.  If  we  substitute  these  three-dimensional  disturbances  into  the 
linearized  three-dimensional  basic  equations,  the  second-order  differential  equation 
for  the  pressure  disturbance,  similar  to  Eq.  2.32,  is  obtained 

2t 


7T  — 


U  —  C 


T  R 
T  +  R 


tt  -  la 2  +  02  -  Ml{u-c)2Cv/R 


RT 


cv/ R  +  (71  -  1) 


>7T  =  0  ,  (4.11) 


where  c  =  u/a  is  the  complex  phase  velocity.  The  following  transformations,  wb'  h 
are  due  to  Squire  [1933],  are  introduced 

=  a2  +  j32  ,  aM\  =  aM\  ,  an  =  an  .  (4.12) 

Thus, 

a  =  a  cos(t?)  ,  0  =  a  sin  (i?)  ,  Mi  =  A/i  cos  (i?)  .  (4.13) 

Using  these  transformations,  Eq.  4.11  can  be  cast  into  the  same  form  as  those  for 
the  flow  with  two-dimensional  wave  disturbances,  t.e., 


n  — 


2  u 

u  —  c 


T  R 
\T  +  R 


n  -  a2  <  1  - 
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RT 
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n  =  0.  (4.14) 
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The  instability  characteristics  of  the  three-dimensional  disturbances 
were  studied  for  compressible  mixing  layers  formed  by  different  perfect  gases.  The 
effect  of  the  angle  of  propagation  direction  d  of  three-dimensional  disturbances  on 
the  instability  behavior  of  mixing  layers  was  considered.  For  the  mean  profiles 
given  by  Eqs.  3.10,  3.15  and  3.20  with  u2/ul  =  0.5,  T2/Ti  =  2.0,  R2/R1  =  0.5 
and  71  =  72  =  1.4,  the  amplification  rate  as  a  function  of  the  propagation  angle 
i9  is  shown  in  Figs.  4.26a  and  4.26b  for  the  cases  Mc\  =  0.25  and  Mc\  =  0.75  re¬ 
spectively.  The  maximum  amplification  rate  —  a,m  decreases  as  the  angle  d  of  the 
disturbance  increases  for  the  low-speed  mixing  layers.  Also,  the  band  of  unstable 
frequencies  becomes  narrower  as  the  angle  1?  of  a  disturbance  increases.  For  the 
high-speed  mixing  layer,  the  maximum  amplification  rate  is  higher  than  that  of  the 
corresponding  two-dimensional  disturbances  over  some  range  of  the  angle  of  prop¬ 
agation  d.  Therefore,  the  three-dimensional  results  may  provide  better  agreement 
with  the  experimental  data  for  high  Mach  number  mixing  layers.  Recall,  however, 
that  according  to  Papamoschou  &:  Roshko’s  experimental  data,  the  growth  rate 
approaches  an  asymptotic  value  as  the  convective  Mach  number  Mc\  becomes  su¬ 
personic.  Even  though  the  growth  rate  of  the  oblique  wave  disturbances  is  higher 
than  the  corresponding  two-dimensional  disturbances,  it  nevertheless  decreases  con¬ 
tinuously  as  the  convective  Mach  number  Mcj  is  increased  (see  Fig.  4.27). 


4.9  Wake-dominated  Mixing  Flows 


Owing  to  the  finite  thickness  of  a  splitter  plate  and  the  displacement 
thickness  of  the  splitter  plate  boundary  layers,  the  initial  region  of  a  two-stream  mix¬ 
ing  flow  is  wake-dominated.  The  mean  quantities  of  the  developing  wake/mixing 
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layer  downstream  of  a  splitter  plate  can  be  obtained  by  numerically  solving  the 
laminar,  two-dimensional,  thin-layer  Navier-Stokes  equations.  In  order  to  simplify 
the  problem,  we  assumed  that  the  mean  velocity  profile  is  composed  of  the  usual 
hyperbolic  tangent  profile  plus  a  wake  component  (due  to  the  splitter  plate)  repre¬ 
sented  by  a  Gaussian  distribution  (see  Eq.  3.12).  For  incompressible  mixing  layers, 
Koochesfahani  &  Frieler  [1987]  studied  the  linear  spatial  instability  characteristics 
of  both  uniform  and  non-uniform  density  plane  mixing  layers,  taking  into  account 
the  wake  component  of  the  initial  velocity  profile.  They  found  that  there  are  two 
unstable  waves  in  incompressible  mixing  layers.  One  is  the  shear  layer  mode  and 
the  other  is  the  wake  mode.  For  compressible  mixing  layers,  the  instability  be¬ 
havior  of  mixing  layers  was  investigated  here  with  a  wake  component  in  the  mean 
velocity  profile.  The  shear  layer  mode  and  the  wake  mode  corresponding  to  that  of 
the  incompressible  mixing  flows  found  by  Koochesfahani  &  Frieler  were  also  found 
for  compressible  mixing  layers.  The  instability  characteristics  of  the  shear  layer 
mode  and  the  wake  mode,  with  different  normalized  wake  deficits  w,  are  shown  in 
Figs.  4.28a  and  4.28b  at  free  stream  Mach  numbers  Mj  =  1  and  Mj  =  2  respectively. 
As  the  normalized  wake  deficit  w  increases,  the  mixing  layer  becomes  more  unstable 
and  the  band  of  amplified  disturbance  frequencies  u>  becomes  narrower  for  both  the 
shear  layer  mode  and  the  wake  mode.  A  comparison  is  made  in  Fig.  4.29  of  the 
normalized  growth  rate  as  a  function  of  the  convective  Mach  number  Mc\  of  both 
shear  layer  mode  and  the  wake  mode  (it;  =  0.4)  with  the  function  F(Mci )  and  with 
the  Papamoschou  &  Roshko  experimental  data.  For  the  shear  layer  mode,  the  nor¬ 
malized  growth  rate  hats  a  similar  behavior  as  the  function  F(Mci ),  but  with  higher 
values  of  the  normalized  growth  rate  in  the  region  0  <  Mcl  <  2.  As  the  convective 
Mach  number  Mc\  increases  continuously  in  the  region  of  supersonic  disturbances, 
the  growth  rate  of  the  shear  layer  mode  starts  to  rise  up  to  a  local  extremum  value 


and  then  falls  again.  For  the  wake  mode,  the  amplification  rate  decreases  very 
slowly  in  the  subsonic  convective  Mach  number  region  and  continuously  decreases 
in  the  region  of  supersonic  convective  Mach  numbers.  Figure  4.30a  shows  the  effect 
of  the  normalized  wake  deficit  w  on  the  normalized  maximum  amplification  rate  of 
the  shear  layer  mode.  If  the  convective  Mach  number  Mc\  is  less  than  about  0.4, 
the  value  of  the  normalized  wake  deficit  w  does  not  affect  the  behavior  of  the  nor¬ 
malized  maximum  amplification  rate  as  a  function  of  the  convective  Mach  number 
Mcl.  But  in  the  region  0.4  <  Mcl  <  1.0,  the  normalized  maximum  amplification 
rate  decreases  faster  for  smaller  values  of  w.  For  supersonic  convective  Mach  num¬ 
bers,  the  value  of  the  convective  Mach  number  Mc\ ,  which  corresponds  to  the  local 
extremum  value,  shifts  to  bigger  values  as  w  decreases.  For  the  wake  mode,  the 
normalized  maximum  amplification  rate  decreases  faster  as  a  function  of  the  con¬ 
vective  Mach  number  Mc\  with  increasing  wake  deficit  w  (see  Fig.  4.30b).  Studies 
of  wake-dominated  mixing  layers  suggest  that  a  larger  value  of  the  growth  rate  at 
large  convective  Mach  number  Mc\  can  be  obtained  by  increasing  the  normalized 
wake  deficit  w. 


4.10  Effect  of  the  Thickness  of  the  Total  Temperature  Profile 


The  thickness  of  the  toted  temperature  profile  may  be  much  smaller 
than  that  of  the  meem  velocity  profile  in  the  initial  region  of  the  flow  near  the  split¬ 
ter  plate  tip.  It  is  interesting,  therefore,  to  find  out  how  this  relative  thickness  can 
affect  the  instability  behavior  of  mixing  layers.  For  incompressible  mixing  layers, 
Koochesfahani  &  Frieler  [1987]  studied  the  effect  of  the  density  profile  thickness  on 
the  instability  behavior  of  mixing  layers.  They  found  that  when  the  high  density 
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fluid  is  carried  on  the  low-speed  side,  the  density  profile  thickness  must  be  smaller 
than  a  certain  value  before  the  wake  mode  becomes  dominant  and  that  when  the 
high  density  fluid  is  on  the  high-speed  side,  regardless  of  the  density  profile  thick¬ 
ness,  the  shear  layer  mode  is  always  dominant.  For  compressible  mixing  layers,  we 
considered  the  effect  of  the  thickness  of  the  total  temperature  profile  relative  to 
the  mean  velocity  profile  on  the  amplification  rate  of  mixing  layers.  In  the  calcula¬ 
tions,  the  hyperbolic  tangent  mean  velocity  profile  given  by  Eq.  3.10  and  the  total 
temperature  profile  given  by  Eq.  3.16  with  ((y)  replaced  by  C(y/<r),  see  Fig.  3.5, 
were  considered,  where  cr  adjusts  the  thickness  of  the  total  temperature  profile  rel¬ 
ative  to  the  mean  velocity  profile.  For  the  velocity  ratio  u2/u]  =  0.5,  the  total 
temperature  ratio  T^/Tn  =  0.5  and  o  =  1.0,  0.8,  0.6,  0.4,  the  variation  of  the 
amplification  rate  —  a*,  as  a  function  of  the  amplified  disturbance  frequencies  u> ,  is 
shown  in  Figs.  4.31a  -  4.31c  for  different  values  of  the  free  stream  Mach  number 
Mi .  As  the  thickness  of  the  total  temperature  decreases,  the  maximum  amplifica¬ 
tion  rate  —  ajm  increases  and  the  band  of  amplified  frequencies  becomes  wider  for 
both  subsonic  and  supersonic  mixing  layers.  The  normalized  maximum  amplifica¬ 
tion  rate  as  a  function  of  the  convective  Mach  number  Mci  is  shown  in  Fig.  4.32  for 
different  values  of  the  relative  thickness  parameter  a.  As  the  thickness  of  the  total 
temperature  profile  a  becomes  smaller,  the  normalized  maximum  amplification  rate 
decreases  slower  in  the  subsonic  convective  Mach  number  region.  For  supersonic 
convective  Mach  numbers,  the  larger  value  of  the  normalized  amplification  rate  is 
obtained  with  decreasing  a. 
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CHAPTER  5 


THE  EFFECT  OF  WALLS  ON  SPATIALLY  GROWING 


SUPERSONIC  MIXING  LAYERS 


The  inviscid  instability,  with  respect  to  supersonic  disturbances,  of  a 
spatially  growing  plane  mixing  layer  inside  parallel  flow  guide  walls  is  investigated 
here  using  linear  stability  analysis.  The  mixing  layer  is  assumed  to  be  inviscid 
and  formed  by  the  same  gases  in  the  two  streams.  The  mean  flow  is  treated  as 
parallel.  The  purpose  of  this  chapter  is  to  give  a  description  of  how  the  instability 
characteristics  of  the  mixing  layer  are  affected  by  the  flow  guide  walls  and  by  the 
distance  between  the  walls.  It  is  found  that  the  maximum  amplification  rates  of  the 
mixing  layers  approach  an  asymptotic  value  and  that  this  maximum  amplification 
rate  increases  to  a  maximum  value  and  decreases  again  as  the  distance  between  the 
parallel  walls  decreases  continuously.  For  a  mixing  layer  inside  parallel  flow  guide 
walls,  the  growth  rate  of  three-dimensional  modes  is  larger  than  the  corresponding 
two-dimensional  mode  at  higher  convective  Mach  numbers.  But  the  growth  rate  of 
two-dimensional  supersonic  instability  waves  is  higher  than  their  three-dimensional 
counterparts  for  a  mixing  layer  inside  a  rectangular  duct  (Tam  &  Hu  [1988,  1989]). 
All  flow  quantities  used  in  this  chapter  are  the  dimensionless  quantities  defined  in 
Chapter  2. 
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5.1  Instability  Behavior  of  Bounded  Mixing  Layers 


The  instability  characteristics  of  mixing  layers  with  mean  velocity 
and  temperature  profiles  given  by  Eq.  3.10  and  Eq.  3.15  are  determined  for  two- 
dimensional,  spatially  growing  disturbances  with  u2ju\  =  0.5,  T2/Ti  =  1.0,  6  =  3 
unit  lengths  (defined  in  Eq.  2.42)  and  d  =  AS  =  12  unit  lengths,  at  different  free 
stream  Mach  numbers  M\.  The  main  result,  shown  in  Fig  5.1,  is  that  when  the 
convective  Mach  number  of  the  flow  is  supersonic,  as  was  also  noted  by  Mack  [1969. 
1984]  in  his  analysis  of  supersonic  boundary  layer  stability,  there  are  many  super¬ 
sonic  unstable  modes,  as  opposed  to  only  one  for  the  case  of  subsonic  convective 
Mach  numbers. 


The  existence  of  many  supersonic  instability  modes  in  the  case  of 
mixing  layers,  has  been  known  from  the  work  of  Tam  &  Hu  [1988,  1989].  Their 
results,  however,  were  based  on  instability  calculations  inside  a  three-dimensional 
rectangular  channel.  Fig.  5.2  (the  normalized  maximum  amplification  rate  vs.  the 
convective  Mach  number  Mc\ )  indicates  that  the  growth  rate  of  the  most  unstable 
supersonic  instability  mode  of  a  mixing  layer  with  two-dimensional,  spatially  grow¬ 
ing  disturbances  approaches  an  asymptotic  value  as  the  convective  Mach  number 
becomes  supersonic,  in  accord  with  the  previously  cited  growth  rate  experiments  by 
Papamoschou  &  Roshko  [1986,  1988].  The  convective  Mach  number  here  is  defined 
by  Eq.  1.4. 


If  a  mixing  layer  inside  parallel  flow  guide  walls  is  subjected  to  a  small. 


51 


three-dimensional,  spatially  growing  disturbance  of  the  form 

q'(x,y,z,t)  =  q(y)  exp  [ i  (  ax  +  (3z  -  w*)],  (5.1) 

where  0  is  a  complex  wave  number  in  the  z  direction  and  a;  is  a  real  disturbance 
frequency,  then  the  growth  rate  of  the  three-dimensional  mode  is  larger  than  the 
corresponding  two-dimensional  mode  at  high  convective  Mach  numbers  over  some 
range  of  propagation  direction  (see  Fig.  5.3).  This  result  is  similar  to  that  of  previous 
studies  of  a  free  mixing  layer  with  three-dimensional  disturbances  by  Gropengiesser 
[1970],  Ragab  &  Wu  [1988],  Jackson  &;  Grosch  [1988]  and  by  Sandham  &  Reynolds 
[1989].  For  a  mixing  layer  inside  a  rectangular  duct,  however,  Tam  &  Hu  [1989] 
considered  three-dimensional,  spatially  growing  disturbances  of  the  form 

P'  =  P(y)  exp  [i  (  ax  -  ut )  ]  cos  (  2mnz/b )  ;  m  =  0,  1,  2,  ...  ,  (5.2) 

where  p'  is  the  perturbation  pressure  and  b  is  the  span  of  the  rectangular  duct.  The 
important  result  of  their  investigations  was  that,  for  reasonably  thick  mixing  layers, 
two-dimensional  supersonic  instability  waves  have  larger  spatial  growth  rates  than 
their  three-dimensional  counterparts.  We  conducted  similar  calculations  as  part  of 
this  investigation.  Our  results  are  in  agreement  with  Tam  &  Hu’s  results.  The 
instability  characteristics  of  three-dimensional  modes  are  characterized  by  similar 
behavior  as  two-dimensional  modes.  As  a  consequence,  two-dimensional  modes  were 
considered  mainly  in  the  work  described  below. 


A  comparison  of  the  current  results  of  growth  rate  with  the  Pa- 
pamoschou  &  Roshko  experimental  results  [1986,  1988]  and  with  the  results  ob¬ 
tained  from  the  two-dimensional  free  mixing  layer  calculations  by  Zhuang  et  al. 
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[1988],  are  made  in  Fig.  5.4.  We  can  see  how  the  parallel  flow  guide  walls  affect  the 
instability  of  the  mixing  layer.  The  existence  of  walls  makes  the  mixing  layer  more 
unstable  and  keeps  the  maximum  amplification  rates  from  reaching  asymptotically 
small  values  for  supersonic  convective  Mach  numbers,  but  has  no  discernible  effect 
on  the  amplification  rate  in  mixing  layers  with  subsonic  convective  Mach  numbers. 


For  a  fixed  high  stream  Mach  number  M\ ,  with  the  same  mean  veloc¬ 
ity  and  temperature  profiles  given  as  described  above,  the  instability  characteristics 
of  mixing  layers  with  two-dimensional  spatially  growing  waves  are  calculated  for 
different  values  of  d  (d  =  20,  16,  12,  8  unit  lengths).  The  most  unstable  super¬ 
sonic  instability  mode,  which  is  not  the  Kelvin-Helmholtz  mode,  moves  to  higher 
frequencies  and  its  maximum  amplification  rate  increases  as  the  distance  between 
two  parallel  walls  decreases  (see  Fig.  5.5).  This  maximum  amplification  rate,  how¬ 
ever,  reaches  its  maximum  value  and  decreases  as  the  distance  between  the  walls 
decreases  continuously  (see  Fig.  5.6).  The  instability  characteristics  of  subsonic 
mixing  layers  are  not  affected  by  the  distance  between  the  walls.  The  thickness  of 
the  mixing  layer  was  chosen  as  the  characteristic  length  in  the  calculations,  which 
made  the  solution  of  the  bounded  mixing  layer,  a s  d/6  —*  oo,  approach  that  of  the 
corresponding  free  mixing  layer.  Note  that  if  the  distance  between  the  parallel  walls 
is  chosen  as  the  characteristic  length,  the  growth  rate  of  the  supersonic  instability 
waves,  scaled  to  this  characteristic  length,  decreases  at  a  fixed  frequency  as  the 
mixing  layer  becomes  thicker  and  thicker  (Tam  &  Hu  [1989]). 
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5.2  The  Pressure  Perturbation  Field 


With  the  computed  eigenfunction  n(y),  the  pressure  perturbation 
fields  can  be  obtained  by 

p'(x,y,t )  =  7r(y)  exp  [  i(ax  -  ut)]  .  (5.3) 

The  real  part  of  p1  represents  the  physical  quantity  of  the  pressure  perturbations. 
Contour  plots  of  the  pressure  perturbation  fields,  which  combine  the  periodic  term 
and  the  growth  term,  are  shown  in  Figs  5.7,  5.8,  5.9  and  5.10.  Note  that  the  con¬ 
vective  Mach  numbers  of  the  mixing  layers  corresponding  to  these  contour  plots 
are  either  close  to,  or  larger  than,  unity,  so  the  growth  of  the  mixing  layers  is 
small  within  the  extent  plotted.  Fig.  5.7  shows  the  flow  with  the  same  supersonic 
convective  Mach  numbers  at  both  sides  of  the  boundaries  (Mc i  =  Mc 2  =  1.375). 
We  can  see  that  the  compression/expansion  waves  propagate  along  the  Mach  an¬ 
gle  vi  =  arcsin(l /Mcl)  or  i/2  =  arcsin(l/Mc2)  for  supersonic  convective  Mach 
numbers.  By  measuring  the  Mach  angle  (u\  =  i/2  =  45.7°),  I  estimated  a  con¬ 
vective  Mach  number  of  A/Cl  =  Mc2  =  1.397,  in  close  agreement  with  Mc\  and 
MC2  estimated  using  Eq.  (1.4).  The  variations  in  the  strength  of  these  compres¬ 
sion/expansion  waves  indicate  that  these  waves  are  reflected  by  the  walls.  Fig.  5.8 
is  the  contour  plot  for  the  free  mixing  layer  with  the  upper  stream  supersonic  and 
the  lower  stream  subsonic  ( Mc\  =  1.81  and  Mc2  =  0.94).  The  upper  flow  compres¬ 
sion/expansion  waves  propagate  with  the  Mach  single  v\  —  34.2°,  i.e.,  Mc  1  =  178, 
with  no  reflections.  From  these  two  contour  plots  we  can  see  that  for  supersonic 
disturbances,  because  of  the  existence  of  the  walls,  the  energy  carried  by  the  wave 
system  is  reflected  back  to  the  mixing  layer  instead  of  being  radiated  and  lost  to  the 
fair  field.  The  feedback  mechanism  between  the  growing  supersonic  mixing  layer  and 
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the  wave  system  makes  the  bounded  supersonic  mixing  layer  more  unstable  than 
the  corresponding  free  supersonic  mixing  layer,  which  is  losing  energy  to  acoustic 
radiation  to  the  far  field  (Mack  [1969],  Tam  [1989]  and  Dimotakis  [1989]).  Figs.  5.9 
and  5.10  Eire  two  other  excimples  of  contour  plots  for  bounded  mixing  layers;  one  is 
with  the  lower  stream  supersonic  and  upper  stream  subsonic,  Eind  the  other  is  vice 
versa.  Reflections  of  the  compression/expansion  wave  system  propagating  along  the 
Mach  angles  for  supersonic  convective  Mach  numbers  can  also  be  seen  in  those  two 
figures  (Zhuang  et  al.  [1989]). 

5.3  Streaklines  of  the  Disturbed  Mixing  Layer 


In  order  to  describe  the  flow  patterns  of  the  mixing  layers,  streaklines 
were  calculated  for  the  modes  corresponding  to  the  maximum  amplification  rates. 
A  streEikline  is  the  locus  of  particles  which  have  earlier  pELSsed  through  a  prescribed 
point.  In  the  experiments,  the  streaklines  can  be  produced  and  visusdized  by  the 
continuous  release  of  marked  particles  (dye,  smoke,  or  bubbles)  from  a  given  point. 
The  procedure  for  the  calculation  of  streaklines  is  the  same  as  that  used  by  Michalke 
[1965].  The  motion  of  a  particle  is  given  by 


dx 

dt 

dy 

dt 


t)u'{x,y,t)  +  u(y)  , 
ej v(x,y,t)  . 


(5.4) 


With  appropriate  initial  conditions,  x(t0)  =  i0  Etnd  y(£0)  =  yo>  &  pathline  can  be 
determined.  In  order  to  plot  each  streakline  for  a  fixed  time,  the  pathlines  for  var¬ 
ious  initial  times,  1 0,  are  calculated.  In  the  calculations,  I  have  chosen  x0  =  0. 
Ci  =  0.0005  for  subsonic  convective  Mach  numbers  and  x0  =  0,  «i  =0.1  for  super- 


sonic  convective  Mach  numbers.  The  results  are  plotted  in  Figs.  5.11  and  5.12  and 
represent  streaklines  corresponding  to  supersonic  and  subsonic  disturbances,  respec¬ 
tively.  These  flow  patterns  show  that,  for  supersonic  disturbances,  the  growth  of 
the  mixing  layer,  as  labeled  by  the  streaklines,  is  much  slower  than  that  for  subsonic 
disturbances,  and  the  spreading  rate  of  the  mixing  layer  is  small. 
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CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 

An  investigation  of  the  linear  instability  behavior  of  compressible  un¬ 
bounded  and  bounded  mixing  layers  with  spatially  growing  wave  disturbances  has 
been  performed.  The  mixing  layer  is  formed  by  two  parallel  streams  with  different 
perfect  gases  and  the  flow  is  assumed  to  be  inviscid  and  non  reacting. 

6.1  Results  for  Unbounded  Mixing  Layers 

For  compressible  unbounded  mixing  layers  with  given  ratios  of  specific 
heats  7i ,  72 ,  the  instability  characteristics  of  flows  were  calculated  in  which  the  free 
stream  Mach  number  Mi  varied  from  0  to  11,  the  velocity  ratio  u^/ui  varied  from 
0.25  to  0.75,  the  temperature  ratio  T2/T1  varied  from  0.2  to  5.0,  the  gas  constant 
ratio  Ri/Ri  varied  from  0.1  to  5.  The  effect  of  compressibility  on  the  properties 
of  mixing  layers  was  determined.  The  sensitivity  to  the  details  of  the  mean  flow 
profiles  on  the  instability  behavior  of  mixing  layers  was  also  studied.  A  summary 
of  the  results  for  unbounded  mixing  layers  is  presented  below: 

1.  For  a  free  mixing  layer  with  hyperbolic  mean  flow  profiles,  there  is 
only  one  unstable  mode  propagating  with  the  phase  velocity  of  the 
most  unstable  wave  C*m  approximately  equal  to  the  isentropic  esti¬ 
mate  of  the  convective  velocity  of  the  large  scale  structures  u*  for 
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subsonic  convective  Mach  numbers.  As  the  convective  Mach  number 
approaches  or  exceeds  unity,  there  are  always  two  unstable  modes. 
One  with  a  phase  velocity  C*m  <  u *  (slow  mode)  and  the  other  with 
a  phase  velocity  C*m  >  u*  (fast  mode). 

2.  For  supersonic  mixing  layers,  the  phase  velocity  of  the  most  unstable 
mode  C;m  decreases  for  the  slow  mode  and  increases  for  the  fast  mode 
with  increasing  free  stream  Mach  number  M\ .  The  bands  of  amplified 
disturbance  frequencies  u  of  these  two  unstable  modes  overlap  for 
some  range  of  frequencies.  For  the  low  supersonic  convective  Mach 
numbers,  the  fast  mode  is  more  unstable  than  the  slow  mode  when 
the  heavy  gas  is  on  the  low-speed  side  and  the  slow  mode  is  dominant 
when  the  heavy  gas  is  on  the  high-speed  side. 

3.  For  a  given  convective  Mach  number,  the  band  of  unstable  frequencies 
u)  becomes  wider  and  the  maximum  amplification  rate  —  atm  becomes 
smaller  as  the  velocity  ratio  u2/ui  increases.  The  phase  velocity  of 
the  most  unstable  mode  C*m  is  a  linear  function  of  the  velocity  ratio 
tt2/ui  for  both  subsonic  and  supersonic  convective  Mach  numbers. 

4.  For  an  incompressible  unbounded  mixing  layer  with  a  hyperbolic  mean 
velocity  profile,  the  maximum  amplification  rate  — ajm  is  approxi¬ 
mately  a  linear  function  of  the  velocity  ratio  A  =  (1  —  u2/ui)/(l  + 
u2/ui)  and  —  a*m  rises  with  increasing  density  ratio  pz/p\  for  all 


0  <  A  <  1  and  density  ratios  P2/P1  <  1.  For  density  ratios  P2/P1 
greater  than  unity,  the  function  — a,m(A)  becomes  nonlinear  and  the 
maximum  amplification  rate  —  a*m  falls  with  increasing  density  ratio 
P2/P1  for  some  values  of  the  velocity  ratio  A. 

5.  The  instability  characteristics  of  unbounded  mixing  layers  are  not  af¬ 
fected  by  the  difference  between  and  72  for  incompressible  flows. 
For  the  subsonic  mode  ( Klc\  7^  0)  and  the  slow  supersonic  mode 
C;m  <  u*,  the  band  of  unstable  frequencies  u>  becomes  wider  and 
the  maximum  amplification  rate  —  atim  rises  if  the  high-speed  stream 
has  the  lower  value  of  the  ratio  of  specific  heats.  But  the  band  of  am¬ 
plified  disturbance  frequencies  becomes  narrower  and  the  maximum 
amplification  rate  falls  for  the  fast  supersonic  mode  C*m  >  u *  when 
the  high-speed  stream  has  the  lower  value  of  the  ratio  of  specific  heats. 

6.  The  eigenfunction  n (y)  decays  exponentially  away  from  the  mixing 
layers  in  the  stream  with  subsonic  disturbances  and  displays  an  oscil¬ 
latory  decay  in  the  stream  with  supersonic  disturbances. 

7.  A  nearly  universal  dependence  of  the  normalized  maximum  amplifica¬ 
tion  rate  on  the  convective  Mach  number  is  found  for  two-dimensional, 
spatially  growing  wave  disturbances,  with  the  normalized  maximum 
amplification  rate  decreasing  significantly  with  increasing  convective 
Mach  numbers  in  the  subsonic  region.  At  high  convective  Mach  num- 
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bers,  the  maximum  amplification  rate  —  o,m  of  mixing  flows  with 
three-dimensional  spatially  growing  disturbances  is  bigger  than  that 
of  corresponding  two-dimensional  disturbances.  But  for  both  two- 
dimensional  and  three-dimensional  disturbances,  the  maximum  am¬ 
plification  rate  decreases  continuously  as  the  convective  Mach  number 
Mc i  increases. 

8.  Adding  a  wake  component  to  the  mean  velocity  profile,  or  decreasing 
the  thickness  of  the  total  temperature  profile  relative  to  the  mean 
velocity  profile,  can  make  the  normalized  maximum  amplification  rate 
decrease  slower  for  both  subsonic  and  supersonic  convective  Mach 
numbers. 

6.2  Results  for  Bounded  Mixing  Layers 

For  compressible  bounded  mixing  layers,  the  inviscid  instability  with 
respect  to  supersonic  disturbances,  of  a  spatially  growing  plane  mixing  layer  inside 
parallel  flow  guide  walls  has  been  investigated  using  linear  stability  analysis.  The 
thickness  of  the  mixing  layer  was  chosen  as  the  characteristic  length  in  the  calcula¬ 
tions,  which  made  the  solutions  of  the  bounded  mixing  layer,  as  d/6  — ►  oo,  approach 
that  of  the  corresponding  unbounded  mixing  layer.  The  effects  of  the  parallel  walls 
and  the  distance  between  the  walls  on  the  instability  characteristics  of  mixing  lay¬ 
ers  were  determined.  Contour  plots  of  the  pressure  perturbation  fields  for  both 
unbounded  and  bounded  mixing  layers  were  calculated.  Also,  the  flow  patterns  of 


the  mixing  layers  were  obtained  by  calculating  the  streaklines  of  the  mixing  layers. 
A  summary  of  the  results  for  bounded  mixing  layers  is  given  below: 

1.  The  existence  of  parallel  flow  guide  walls  makes  the  mixing  layer 
more  unstable  and  keeps  the  maximum  amplification  rate  —  a,m  from 
reaching  an  asymptotically  small  value  for  supersonic  convective  Mach 
numbers,  but  has  no  discernible  effect  on  mixing  layers  with  subsonic 
convective  Mach  numbers. 

2.  The  most  unstable  supersonic  instability  mode,  which  is  not  the 
Kelvin-Helmholtz  mode,  moves  to  higher  frequencies  and  its  maxi¬ 
mum  amplification  rate  — aim  increases  to  its  maximum  value  and 
decreases  again  as  the  distance  between  the  walls  decreases  continu¬ 
ously. 

3.  For  supersonic  convective  Mach  numbers,  the  reflections  of  the  com¬ 
pression/expansion  waves  caused  by  the  parallel  walls  provide  a  feed¬ 
back  mechanism  between  the  growing  mixing  layer  structures  and  the 
wave  system.  Bounded  mixing  layers  axe  more  unstable  than  the 
corresponding  unbounded  mixing  layers.  The  maximum  simplifica¬ 
tion  rates  —  atjm  of  the  bounded  mixing  layers  are  found  to  approach 
an  asymptotic  value  for  supersonic  convective  Mach  numbers.  The 
streaklines  of  the  flow  confirm  that  the  spreading  rate  of  the  mixing 
layer  is  unusually  small  for  supersonic  disturbances. 
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Fig.  3.1 


Hyperbolic  tangent  mean  velocity  profiles  for  different  values  of  the  veloc¬ 
ity  ratio  ui/ui. 
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Fig.  3.2 


Mean  velocity  profiles  for  different  values  of  the  normalized  wake  deficit 
w  at  the  velocity  ratio  U2/U1  =  0.5. 
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Fig.  3.3  Mean  temperature  profiles  for  different  values  of  the  temperature  ratio 
T2/T1  for  the  case  u2/ui  =  0.5,  R2/R\  =  1.0  and  M\  =  5.0. 
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Fig.  3.4  Mean  temperature  profiles  for  different  values  of  the  free  stream  Mach 
number  M-y  for  the  case  u^/xiy  =  0.5,  T^/Ti  =  1.0  and  i?2/-Ri  =  1.0. 
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Fig.  3.5  Hyperbolic  tangent  total  temperature  profiles  for  different  values  of  the 
interface  thickness  a  at  the  total  temperature  ratio  Ta/Tn  =  0.5. 


74 


Fig.  3.6  Hyperbolic  tangent  mean  temperature  profiles  for  different  values  of  the 
temperature  ratio  T2/T1. 


Fig.  3.7  Mean  gas  constant  profiles  for  different  values  of  the  gas  constant  ratio 

R2/R1. 
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Fig.  4.4  Phase  velocity  of  the  most  unstable  mode  Cpm  and  the  convective  velocity 
uc  vs.  the  free  stream  Mach  number  M\  for  subsonic  disturbances. 


Fig.  4.5  Instability  characteristics  for  the  case  U2/U1  =  0.25,  T^/Ti  =0.5,  R2/R1  = 
1.0  and  Mi  =  2.25  at  supersonic  disturbances,  (a)  The  amplification  rate 
—a,,  (b)  The  phase  velocity  Cp. 
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Fig.  4.7  Instability  characteristics  for  the  case  U2/U1  =0.25,  T2/T1  =  1.5,  R2/R1  = 
1.0  and  M\  =  3  at  supersonic  disturbances,  (a)  The  amplification  rate 
— a*,  (b)  The  phase  velocity  Cp. 
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Fig.  4.9  Two-dimensional  free  mixing  layer  with  supersonic  disturbances:  U2/U1  =  # 

0.5,  T2/T1  =  1.0  and  R2/R1  =  5.0.  (a)  The  phase  velocity  of  the  most 
unstable  mode  Cpm  and  the  convective  velocity  uc  vs.  the  free  stream 
Mach  number  Mi.  (b)  The  maximum  amplification  rate  — Oim  vs.  the 
free  stream  Mach  number  Mi . 


Fig.  4.10  Two-dimensional  free  mixing  layer  with  supersonic  disturbances:  u^/ux  — 
0.25,  T2/T1  =  1.5,  R2/R1  =  10.  (a)  The  phase  velocity  of  the  most 
unstable  mode  Cpm  and  the  convective  velocity  uc  vs.  the  free  stream 
Mach  number  M\.  (b)  The  maximum  amplification  rate  — o,m  vs.  the 
free  stream  Mach  number  M\. 
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Fig.  4.11  The  maximum  amplification  rate  — a,m  vs.  X  =  (1  —  ui/ui)l(\ -I- u2/«i ) 
for  the  case  T2/T1  =  1.0  and  R2/R1  =  1.0. 
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Fig.  4.12  The  amplification  rate:  T2/T1  =  0.5  and  R.2/R\  =  1.0.  (a)  The  sub¬ 
sonic  mode  Mci  =  0.  (b)  The  subsonic  mode  Mci  =  0.3.  (c)  The  slow 
supersonic  mode  Mci  =  1-36.  (d)  The  fast  supersonic  mode  Mc 2  =  1.51. 
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Fig.  4.13  Phase  velocity  of  the  most  unstable  mode  Cpm  and  the  convective  velocity 
uc  vs.  the  velocity  ratio  u2/iZ]  for  the  case  T2/Ti  =  0.5  and  R2/R1  =  1.0. 
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FIG.  4.14  Growth  rate  plotted  against  the  density  ratio  pi/*,  for  the  case  u^Ul 

0. 
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Fig.  4.16  Growth  rate  plotted  against  A  =  (1  -  u2/t/i)/(l  4-  ui/u\). 


Fig.  4.18  The  amplification  rate:  u^/u^  =  0.5,  T2/T1  =  1.0  and  R2/R1 
(a)  The  subsonic  mode  Mc  1  =  0.0.  (b)  The  subsonic  mode  A/Cl  : 


Fig.  4.21  Eigenfunction  of  the  pressure  disturbances  n(y)  =  nr(y)  +  iVj(y)  for 
supersonic  disturbances:  u2/tti  =  0.25,  T2/Ti  =  0.5  and  R2/R1  =  1.0. 
(a)  The  supersonic  slow  mode  Mi  =  2.6.  (b)  The  supersonic  fast  mode 
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Fig.  4.22  Continued. 
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FlG.  4.23  Least  squares  fitting  the  normalized  maximum  amplification  rate  vs.  the  convective  Mach  number  Mci. 
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Fig.  4.24  Normalized  maximum  amplification  rate  vs.  Mc\  for  hyperbolic  tangent 
mean  temperature  profiles  comparison  with  F{MC\ ). 
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Fig.  4.25  A  comparison  of  F{MC\ )  with  Ragab  &  Wu’s  numerical  data  and  with 
Papamoschou  &  Roshko’s  experimental  data. 
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Fig.  4.27  Normalized  growth  rate  vs.  Mc\. 
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FIG.  4.28  Instability  characteristics  of  the  shear  layer  mode  and  the  wake  mode  for 
different  values  of  the  normalized  wake  deficit  w.  u2/“i  =  0.5,  Ti/T\  = 
1.0  and  R2/Ri  =  1.0.  (a)  M,  =  1.0.  (b)  Mx  =  2.0. 


Fig.  4.30  Normalized  maximum  amplification  rate  vs.  Mci:  U2 /u\  —  0.5,  Tj/Tj 
1.0  and  R2/R1  =  1.0.  (a)  The  shear  layer  mode,  (b)  The  wake  mode. 
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Fig.  4.31  Effect  of  the  thickness  of  the  total  temperature  profile  a  on  the  am¬ 
plification  rate  of  the  mixing  layers:  u2/ui  —  0.5,  Tt2/Tt i  =  0.5  and 
R2/R1  =  1.0.  (a)  Mi  =  0.0.  (b)  Mi  =  4.0.  (c)  Mx  =  6.0. 
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FIG.  4.32  Normalized  maximum  amplification  rate  ua.  Mci  for  the  case  U2/U1  = 
Tt2/Tt  1  =  0.5  and  R2/R1  =  1.0. 
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Fig.  5.2  Normalized  maximum  amplification  rate  vs.  Mc  1  at  d  =  12  unit  lengths 
and  6  =  3  unit  lengths. 


Fig.  5.3  Instability  characteristics  for  3-D  spatially  growing  disturbances  for  the 
case  u2/ui  —  0-5,  T2/T1  =1.0  and  M\  =  5.5  at  d  =  12  unit  lengths  and 
6  =  3  unit  lengths. 
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Fig.  5.4  A  comparison  of  the  data  for  2-D  bounded  mixing  layer  with  Papamoschou 
&  Roshko’s  experimental  data  and  with  the  function  F(Mci),  which  was 
obtained  by  least  squares  fitting  the  data  from  the  calculations  of  2-D 
mixing  layers. 
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ressure  perturbation  field  of  the  most  unstable  mode  for  the  free  mixing 


Streaklines  for  the  case  u2/«i  =  0.5  and  T2/Ty  =  1.0  at  6  =  3  unit  lengths,  d  =  12  unit  lengths  and 
Mi  =  5.5  (supersonic  convective  Mach  numbers  Mc\  =  Mc2  =  1.375). 
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Fig.  5.12  Streaklines  for  the  case  U2/U1  =  0.5  and  Ti/T\  =  1.0  at  6  =  3  unit 
lengths,  d  =  12  unit  lengths  and  M\  =  2.0  (subsonic  convective  Mach 
numbers  Mcl  =  Mc2  =  0.5). 
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